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ABSTRACT 

This is the fourth paper in a series that reports on our investigation of the clustering properties 
of active galactic nuclei (AGNs) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital 
Sky Survey (SDSS). In this paper we investigate the cause of the X-ray luminosity dependence of the 
clustering of broad-line, luminous AGNs at 0.16 < z < 0.36. We fit the Ha line profile in the SDSS 
spectra for all X-ray and optically selected broad-line AGNs, determine the mass of the supermassive 
black hole (SMBH), Mbh, and infer the accretion rate relative to Eddington (L/Ledd). Since Mbh 
and L/Ledd are correlated, we create AGN subsamples in one parameter while maintaining the same 
distribution in the other parameter. In both the X-ray and optically selected AGN samples, we detect 
a weak clustering dependence with Mbh and no statistically significant dependence on L/Ledd- We 
find a difference of up to 2.7cr when comparing the objects that belong to the 30% least and 30% most 
massive A/bh subsamples, in that luminous broad-line AGNs with more massive black holes reside 
in more massive parent dark matter halos at these redshifts. These results provide evidence that 
higher accretion rates in AGNs do not necessarily require dense galaxy environments, in which more 
galaxy mergers and interactions are expected to channel large amounts of gas onto the SMBH. We 
also present semianalytic models that predict a positive Mdmh dependence on Mbh, which is most 
prominent at Mbh ~ 10®“® Mq. 

Subject headings: galaxies: active - cosmology: large-scale structure of universe - X-rays: galaxies 


1. INTRODUCTION 

There has been increasing interest in large-scale clus¬ 
tering measurements of active galactic nucleus (AGNs) 

in recent year s. _ Such measurements (see review by 

iKrumoe et al.l 1201411 not only allow one to study the 
distribution of matter in the Universe out to redshifts 
where it becomes very challenging and observation- 
ally expensive to detect galaxies in large numbers, but 
they also can be used to constrain theoretical mod¬ 
els of AGN/galaxy coevolution, feedback mechanisms, 
AGN host-galaxy properties, the distribution of AGNs 
with dark matter halo (DMH) mass, and the fueling 
process(es) of supermassive black ho l es (SMBH) (e.g. 
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of thousands of galaxies yield significant clustering de¬ 
pendences on galaxy properties such as luminosity, 
morphological type , spectral type, and redshift (e.g., 
Norberg et aLll2002l: IMadgwick et ^120031; IZehavi et alJ 

200561120llHMenenx et al.ll2006ll2nn9HCoil et al.llfM 7 

These findings confirm the hierarchical model of struc¬ 
ture formation in which more massive, and hence more 
luminous, galaxies reside in more massiv e DMHs and 
are th e refore clustered more st rongly (e.g.. IZehavi et aP 
l2005bL 120 lit ICoil et all 1200611 . Whether this relation 
should also apply to AGN luminosity is not trivial. Since 
the AGN luminosity depends primarily on the mass of 
the SMBH (Mbh), the accretion rate relative to Edding¬ 
ton (L/Ledd), and the radiative accretion efficiency, a 
relation between clustering and AGN luminosity should 
ultimately be due to a connection between DMH mass 
and one or more of these physical parameters. 

Based on sm oothed-particle h ydrodynamic simu l ations 
using GADGET (iSpringell 1200511 . iBooth fc Schavd (1201611 
explore the correlation between SMBH mass and the 
mass of the hosting DMH. In the simulations, the black 
holes grow either by the accretion of ambient gas or 
mergers. The black holes also inject a fixed fraction 
of the rest mass energy of the gas into the surround¬ 
ing medium. A self-regulating black hole injects enough 
energy to displace gas from the host galaxy on longer 
timescales. The binding energy of the gas is dete r mined 
by the DMH potential. Thus, IBooth fc Schavel ()20in[ l 
conclude that the mass of the SMBH is regulated pri¬ 
marily by t he DMH mass and no t the stellar mass of 
the galaxy. iFanidakis et al.l (1201211 us e the semianalyt- 
ical galaxy-formation model GALFORM (iGole et al.l 120001 : 
more details on their simulations are given below in 
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Sect. l6.4|) and find a correlation between SMBH mass and 
DMH mass at almost all cosmic times (z = 0.0 — 6.2). 
iVolonteri et al.l (|2011ll use the measured black hole mass, 
velocity dispersion cr, a nd asymptotic cir c ular v elocity of 
25 local galaxies from iKormendv et al.l (1201111 to show 
that, although with some scatter, the black hole masses 
correlate well with the parent DMH masses. 

AGN clustering studies can be used to observationally 
test such predictions. Several studies have measured the 
large-scale clustering dependence on AGN on properties 
such a s X-ray luminosity, with varying results. IGoil et SI 
(|2009ll find no correlation between the AGN X-ray lumi- 
nosi ty and the clustering strength at z = 0.7—1.4. At z ^ 
0.1, iMountrichas fc Georgakakli (|2012ll also detect no 
significant correl ation in thei r samp le of XMM-Newton- 
selected AGNs. lYang et all (I2006ll reports a tentative 
('^Icr) dependence of the clustering signal, such that 
the brig hter sample is more clu stere d than the fainter 
sampl e. iGappelluti et al.l (1201011 and iKoutoulidis et al.l 
(|2013ll verify this fi nding using differe n t AGN samples at 
different redshifts. IKoutoulidis et al.l (1201311 use a large 
sample of ^1500 AGNs from the Ghandra Deep Field 
(CDF) North and South, the extended CDF South, COS¬ 
MOS, and AEGIS. The median redshifts of their low and 
high Lx AGN samples are (z) ~ 0.8 — 1.1. Except for 
the CDF South survey, they find weak X-ray luminosity 
dependences at a level of up to ^2cr. 

One way to reduce the uncertainties involved in cluster¬ 
ing measurements — beyond using larger AGN samples 
— is to compute the cross-correlation function (CCF) 
with a dense galaxy sample instead of computing the 
autocorrelation f unction (ACF) of the AGNs. Sev¬ 
eral studies (e.g. , iLi et all [Mia IGoil et al.l[20nl IMil: 
Wake et al.ll200^ iHickox et al.ll2009l IMountrichas et al l 
2013HGeorgakakis et al.ll2014ll demonstrate the potential 
of this approach and compute the CCF between AGNs 
and a large sample of galaxies to infer the ACF of the 
AGNs. The significant increase in the number of pairs 
at a given separation, used to measure the clustering 
strength, reduces the uncertainties in the spatial corre¬ 
lation function compared to the direct measurements of 
the AGN ACF. 

In iKrumpe et al.l (l2010al hereafter P aper 1), we use 
the same technique as IGoil et al.l (1200911 to measure the 
CCF between ROSAT All-Sky Survey (RASS) AGNs 
identified in the Sloan Digital Sky Survey (SDSS) and 
a large set of SDSS luminous red galaxies (LRGs) at 
0.16 < z < 0.36. The study is based on SDSS data 
release 4-1- (DR4-I-). The unprecedented low uncertain¬ 
ties of the inferred broad-line AGN ACF allow us to split 
the sample into low and high X-ray luminosity subsam¬ 
ples and to report a '^2.5cr X-ray luminosity dependence 
of broad-line AGN clustering. We find that higher lumi¬ 
nosity AGNs cluster more strongly than their lower lu¬ 
minosity counterparts. Consequently, higher luminosity 
X-ray AGNs reside, on average, in more massive DMHs 
than do lower luminosity X-ray AGN^ __ 

In the second paper of this series (IMivaii et al.ll20lTL 
hereafter Paper II), we describe a novel method of ap¬ 
plying the halo occupation distribution (HOD) modeling 
technique directly to the precise measured CCF between 
RASS/SDSS AGNs and SDSS LRGs to constrain the dis¬ 
tribution of AGNs as a function of DMH mass. This 
method also allows us to derive the large-scale bias pa¬ 


rameter of the AGN sample with much lower systematic 
uncertainties than using a phenomenological power-law 
fit, as is often done. As shown in Paper H, the X-ray lu¬ 
minosity dependence is more prominent in the one-halo 
term (rp < 1 h~^ Mpc). The HOD-based typical DMH 
masses derived from the two-halo term for the high- and 
low-luminosity RASS/SDSS AGN subsamples differ by 
^l.StT. In addition, we find that models where the AGN 
fraction among satellites decreases with DMH mass be¬ 
yond Mdmh 10^^ h~^ Mq are preferred. This is in 
contras t to what is found for satellite galax ies without 
AGNs (jZheng et al.ll2009l: iZehavi et al.ll2011[) . 

In the third paper ( Krumpe et al.ll20I2L hereafter Pa¬ 
per HI), we extend the cross-correlation measurements to 
lower and higher redshifts, covering a redshift range of 
z = 0.07 — 0.50. We show that the weak X-ray luminos¬ 
ity dependence of broad-line AGN clustering is also found 
if radio-detected AGNs are excluded. Furthermore, we 
compute the large-scale clustering for optically selected 
broad-line SDSS AGNs using the final SDSS DR7, but we 
detect no optical luminosity dependence on the cluster¬ 
ing strength, although the optical broad-line SDSS AGN 
sample in 0.16 < z < 0.36 contains more than twice as 
many objects as the RASS/SDSS AGN sample. We con¬ 
clude that the most likely explanation for this result is 
the smaller dynamic range probed in optical luminosities 
compared to X-ray luminosities. 

In this paper we focus on the redshift range 0.16 < z < 
0.36, where the CCF of RASS/SDSS AGNs and LRGs 
has the highest signal-to-noise ratio (S/N). We fit the 
Ha line profile in the SDSS spectra of broad-line AGN 
to infer the A/bh and the L/Ledd- Dividing the AGN 
sample into low and high Mbh mass, as well as low and 
high L/Ledd, provides insights into the main physical 
driver of the weak detected X-ray luminosity dependence 
of broad-line AGN clustering. 

This paper is organized as follows. In Section 2, we 
describe the properties of the LRG tracer set and the 
AGN samples. Section 3 provides details on how we 
fit the Ha line profile in the optical SDSS AGN spec¬ 
tra, derive the Mbh, estimate L/Ledd, and define our 
AGN subsamples. In Section 4, we briefly summarize 
the cross-correlation technique, how the AGN ACF is in¬ 
ferred from this, and how we derive the large-scale bias 
parameters using HOD modeling. Section 5 provides the 
results of our clustering measurements. The detailed re¬ 
sults of the HOD modeling of the CCFs presented in this 
paper and in paper HI will be included in a future pa¬ 
per (T. Miyaji et al. in preparation). Our results are 
discussed in Section 6, and we present our conclusions 
in Section 7. Throughout the paper, all distances are 
measured in comoving coordinates and given in units 
of h“^Mpc, where h = i7o/100kms“^ Mpc“^, unless 
otherwise stated. We use a cosmology of Dm = 0.3, 
Da = 0.7, and as{z = 0) = 0.8, whi ch is consistent wit h 
the WMAP data, release 7 (Table 3 of iLarson et al.l[20llll . 
The same cosmology is used in Papers I-III. Luminosi¬ 
ties and absolute magnitudes are calculated for h = 0.7. 
We use AB magnitudes throughout the paper. All uncer¬ 
tainties represent Icr (68.3%) confidence intervals unless 
otherwise stated. 


2. DATA 
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The data sets used in this study are drawn from the 
SDSS, which consists of an imaging survey in five bands 
and an extensive spectroscopic follow-up survey with a 
fiber spectrograph. The selection of the optically s electe d 
AGN candidates is described in [Richards et al.l (1200211 . 
LRGs are chosen by following lEisenstein et al.l (120011) . 

2.1. SDSS Luminous Red Galaxy Sample 

The selection of SDSS LRGs follows the procedure de¬ 
scribed in Section 2.1 of Paper I and Section 2.2 of Paper 
III. Here we briefly summarize the sample selection. We 
extract LRGs from the web-based SDSS Catalog Archive 
Server Jobs SysteirQ using the flag “gala xy_red.” which is 
based on the selection criteria defined in lEisenstein et al] 
(j2001fl . We verify that the extracted objects meet all 
LRG selection criteria and create a volume-limited sam¬ 
ple with 0.16 < z < 0.36 and —23.2 < < —21.2, 

where is based on the extinction-corrected r*gti:o 

magnitude, ^-corrected and passively evolved to rest- 
frame 5petro niagnitudes at z = 0.3. We consider only 
LRGs that fall into the SDSS area with a DR7 spectro¬ 
scopic completeness ratio of greater than 0.8 and that 
have a redshift confidence level of greater than 0.95. The 
SDSS geometry and compl eteness ratio are expressed i n 
terms of spherical polygons (jHamilton fc Tegmarld[2003l . 
The file is publicly availabl430 through the New York 
Universi ty Value-Added Gal axy Gatalog (NYU-VAGG) 
website (jBlanton et al.ll^OOSf) . 

We co rrect for the SPSS fiber colli sion as described in 
detail in lKrumne et al.l (j2010aL [201211 . We have to assign 
to approximately 2% of all LRGs in our sample a redshift 
due to the fiber collision problem. 

The construction of the random LRG sample is iden¬ 
tical to the procedure described in Section 3.1 of Pa¬ 
per 1. Our LRG random sample contains 200 times as 
many objects as the real LRG sample. We generate a 
set of random R.A. and deck values within DR7 areas 
with spectroscopic completeness ratios of greater than 
0.8, populate areas with higher completeness ratios more 
than ones with lower completeness ratios, and randomly 
assign redshifts to the objects in the sample by using the 
smoothed redshift profile of the observed redshift distri¬ 
bution. 

2.2. RASS/SDSS AGN Samples 

The ROSAT All-Sky Survey iRASS. IVoges et aI1ll999H 
is currently still the most sensitive al l-sky survey in 
the soft (0. 1-2.4 keV) X-ray regime. lAnderson et al.l 
(|2003l 1200711 positionally cross-correlate RASS sources 
with SDSS spectroscopic objects and classify RASS- and 
SDSS-detected AGNs based on SDSS DR5. They find 
6224 AGN with broad permitted emission lines in ex¬ 
cess of 1000 kms“^ FWHM and 515 narrow permitted 
emission line AGNs matching RASS sources within 1 ar- 
cmin. More details on the sample selection are given i n 
Section 2.2 of Paper I and lAnderson et ^ (120031 . |200^ . 
Since ROSAT observed the sky in the soft energy band 
(0.1-2.4 keV), the RASS/SDSS AGN sample is biased to¬ 
ward AGNs with little to no X-ray absorption. The vast 

^ http://casjobs.sdss.org/CasJobs/ 

® http://sdss.physics.nyu.edu/lss/dr72 
® http://sdss.physics.nyu.edu/lss/dr4plus 


majority of the optical counterparts are therefore AGNs 
with broad emission lines and UV excess. There is no 
overlap between the RASS/SDSS AGNs and the LRG 
sample. 

To study the X-ray luminosity dependence of the clus¬ 
tering of these AGNs, we have to limit the SDSS foot¬ 
print to the publicly available DR4-I- geometry. We split 
our sample in the redshift range 0.16 < z < 0.36 into 
subsamples according to X-ray luminosity. We use a 
0.1-2.4 keV observed luminosity cut (assuming a pho¬ 
ton index P = 2.5, corrected for Galactic absorption) 
of log (Lx/[ergs“^]) = 44.29. The observed flux in 
the 0.1-2.4 keV band has a large soft-excess contribu¬ 
tion that is not representative of the underlying intrin¬ 
sic, hard power-law X-ray spectrum. Thus, we use the 
template XMM-New ton spectrum of powerfu l radio-quiet 
ROSAT QSOs from iKrumoe et all (l2010bf) to estimate 
the corresponding flux in the 0.5-10 keV and 2-10 keV 
energy ranges. Using the median redshift of the sam¬ 
ple (z = 0.27), the cut at 0.1-2.4 keV corresponds 
to a cut at log (Lo.5-iokev/[ergs“^]) = 43.7 and log 
(L2-iokev/[ergs“i]) = 43.4. 

As an alternative estimate, we mat ch our RASS AGN 
sample with the 3XMM-DR5 catalog ([Rosen et al.ll201^ . 
We fit a regression line between the RASS 0.1-2.4 keV 
fluxes and the XMM-N ewton 2-12 keV fluxes. Using 
Xspec ([Arnaud al.l 119961 1. we model the same ratio with 
a broken power law and determine the luminosity ratios 
in the ranges 0.1-2.4 keV, 0.5-10 keV, and 2-10 keV. We 
find somewhat smaller corrections: log (Lx/[ergs“^]) = 
44.29 corresponds to log (Lo. 5 -iokev/[ergs“^]) = 43.9 
and log (L 2 -iokev/[ergs“^]) = 43.7. This approach has 
the disadvantage that the observations were taken over 
a decade apart, and temporal variation in the X-ray 
luminosity/flux of the objects will affect the estimate. 
However, one could argue that such variations are effec¬ 
tively averaged over in a large sample of X-ray objects 
as used here. Fits to the 0.5-10 keV XMM-Newton spec¬ 
tra verified that the vast majority (^95%) of the cross- 
matched RASS/XMM AGNs are unabsorbed in the X- 
rays (A^h < 10^^ cm“^), while the remaining sources 
show absorption at a level of only a few 10^^ cm“^. 

We calculate the comoving number densities as de¬ 
scribed in detail in Paper 1. For a given R.A. and deck, we 
compute the limiting observable RASS count rate and in¬ 
fer the absorption-corrected flux limit versus survey area 
for RASS/SDSS AGNs. We then compute the comoving 
volume avai lable to each object ( 14) to be included in 
the sample (jAvni fc Bahcalll 1 198011 . The comoving num¬ 
ber density follows by computing the sum of the available 
volume over each object, ttagn = 1/^,*- The comov¬ 

ing number densities for the total, high Lx, and low Lx 
subsamples are 6.0, 0.12, and 5.8 xl0“^ Mpc“^, re¬ 
spectively. 

3. DERIVING Mbh AND L/Ledd THROUGH Ha LINE 
PROFILE FITS 

The so-called virial method is routinely employed to 
estimate Mbh from single-epoch spectra based solely 
on the width of broad emission li nes together with 
an appropriate cont i nuum l uminosity (iKasoi et akll2000l : 
iPeterson fc Wand^ 120001: iVestergaardl 1200^ The 
FWHM of H/3 (FWHMh/ 3 ) and the rest-frame contin- 
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uum luminosity at 5100 A (L 5100 ) are most commonly 
used to compute Mbh for low-redshift AGNs, because 
the corresponding calibrations are the best-studied ones 
to date. On the other hand, the Ha line offers a higher 
S/N than H/3 and tight relations between FWHMh /3 
and FWHMna, as w ell as between L.r,inn and Lurv, 
have been established (|Greene fc Holl200a . iSchulze et al.l 
IMoli . Consequently, accurate Mbh estimates can also 
be obtained from the broad Ha line. 

The simultaneous use of the broad Ha line as a proxy 
for the velocity dispersion in the broad-line region (BLR) 
and for the AGN luminosity has two major advantages. 
First, we avoid any potential contamination of L5100 by 
host galaxy light and the influence of the broad Fe II 
emission i n type I AGNs on th e broad H/3 FWHM mea¬ 
surement (jGreene fc Hd[2n?)^ . Second, we obtain reli¬ 
able BH mass estimates even for those objects where the 
S/N of H/3 would be too low. 

3.1. Fitting the Ha Line Profile 

We retrieve the individual SDSS DR7 optical spectra 
for our entire RASS/SDSS AGN sample in order to mea¬ 
sure the luminosity and FWHM of Ha for each AGN. 
Although the FWHM and integrated flux of Ha could in 
principle be measured directly from the spectra, the nar¬ 
row Ha and [N II] lines on top of the broad Ha line 
have a significant impact on the estimated Mbh and 
L/LHnn valu e s (Fig . [U left panels), as highlighted by 
iDennev et al.l (|2009ll . The narrow Ha originates from 
star formation in the host galaxy and the narrow-line re¬ 
gion (NLR) photoionized by the AGN outside of the cen¬ 
tral. Deblending the broad and narrow emission lines is 
therefore essential to minimize systematic effects by sepa¬ 
rating the contribution from the BLR and the NLR/host 
in the Ha line. 

One major complication in modeling the AGN spec¬ 
tra is the commonly complex, non-Gaussian shape of 
the broad Ha line. Previous s tudies have used multi¬ 
ple Gaussian c ompo n ents fe.g.. iBrotherton et al.l[T994 
Sulentic et al.l 120021 : iShen et al.l 120081 : LSchulze et al.1 
200l~ or high-order Gauss-Hermite polynomials (e.g., 
Salyiander et al.l [20071 : iMcGill et al.ll200^ iStern fc Laorl 
2 OI 2 II to describe the asymmetries of the line. We test 
both approaches to fit the SDSS spectra and conclude 
that the latter is more robust against the choice of initial 
parameter estimates. 

We simultaneously model the broad H/3 and Ha lines 
together with the narrow H/3, [O III] A4960, A5007, Ha, 
and [N II] A6548, A6583 emission lines. We briefly de¬ 
scribe our algorithm of modeling th e broad Ha line below 
and follow an approach similar to IStern fc Laorl (j20I2f ) . 
The model of the H/3 region (left panels in Fig. [T]) is 
only used to better constrain the parameters of narrow 
lines. The well-isolated narrow [O III] line in the H/3 re¬ 
gion gives parameter restrictions that are superior to the 
narrow-line profiles in the Ha region, where the narrow 
and broad lines can be strongly blended. 

As a first step, the approximate underlying continua 
are independently subtracted from the H/3 and Ha lines 
by linear interpolation between the adjacent rest-frame 
spectral regions 4720-4760 A and 5070-5130 A, as well 
as 6150-6230 A and 6750-6950 A, respectively. The fit 
requires initial guesses for the broad and narrow line 


SDSS J102812.30+660723.9 (z=0.198) 



SDSS J172711.80+632241.8 (z=0.217) 



SDSS J132352.75-020101.7 (z= 0.223) 



SDSS J004319.74+005115.4 (z=0.308) 




Fig. 1.— Observed SDSS spectra (black solid lines) and the corre¬ 
sponding models (red solid lines) in the observed H/3 (left panels) 
and Ho (right panels) wavelength ranges. The fitted individual 
line components are indicated by blue solid lines (underlying con¬ 
tinuum, narrow and broad lines). The most prominent line features 
are labeled. 


components. The initial guesses for the FWHM and 
line flux of the broad lines are determined by a fourth- 
order Gauss-Hermite polynomial of the entire Ha line 
and assuming standard Ha and H/3 ratios. The initial 
guesses of the narrow line components are based on the 
total flux of the observed [O III] line and typical narrow- 
line AGN ratios. An eighth-order Gauss-Hermite poly- 
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nomial provides a good fit for the full model of com¬ 
plex broad H/3 and Ha line shapes in almost all cases. 
All narrow emission lines are assumed to have ordinary 
Gaussian profiles that are coupled in redshift, velocity 
dispersion, and line ratios for doublets. While [O III] 
can e xhibit a blueshifted w ing in cases with an outflow 
(e.g.. iMullanev 0111111201311 . the narrow component usu¬ 
ally strongly dominates the line flux, and therefore a 
single-component line fit is generally robust when cou¬ 
pled to lower ionization lines like [N II]. In cases where 
the outflow component dominates, we decouple [O III] 
from the fit to other narrow lines, as described below. 
We also add the two dominant Fe II lines in the 
wavelength range as Gaussians to our model. The best- 
fit model is determined by using a Levenberg-Marquardt 
minimization algorithm, where spectral regions of uncon¬ 
sidered weak narrow emission lines, e.g., [O I]A6300 and 
[S II]A6718, A6732, have been masked out along with pix¬ 
els assigned as bad by the SDSS spectroscopic pipeline. 

We deviate slightly from the scheme above in the fol¬ 
lowing three cases: (1) the broad H/3 line is too weak to 
be modeled, (2) the [O III] line exceeds a width of 600 
km sec“^, or (3) the narrow Ha line flux is negative in 
the initial best-fit model. In the first case, we repeat 
the model without any broad H/3 and Fe II components. 
In the second and third case, the [O III] line is likely to 
be affected by outflows and may not match the line pro¬ 
file of the other narrow lines. Therefore, we allow the 
[O III] line to have an independent redshift and velocity 
dispersion with respect to the other narrow lines. 

One percent of the RASS/SDSS spectra do not allow 
such an analysis because significant parts of the spectral 
range around Ha or [O III] are masked out. We do not 
expect any significant impact on our study, due to the 
low number of objects for which this is the case (see also 
Sect. [5T] below). 

Examples of the best-fit models for typical S/N AGN 
spectra are shown in Figure [T] From those best-fit mod¬ 
els we measure the FWHM and integrated flux of the 
broad-line Ha component, as these are the parameters 
required to estimate Mbh and L/Ledd- Our model is 
not designed to provide a perfect fit to Ha line profiles 
for the most complex cases. However, the model devia¬ 
tions concerning the Ha width and flux are much smaller 
than t he systematics of viri al Mbh estimates (see discus¬ 
sion in iDennev et al.l l2009f) . To assess the uncertainties 
of the derived parameters and to identify unreliable mod¬ 
els, we generate 100 realizations of the same spectrum by 
replacing each pixel within its Icr variance given by the 
SDSS error spectrum and analyze the artificial spectra 
in the same manner. The uncertainty for each quantity 
is then taken as the dispersion (Icr error corresponding 
to 68.3%) in the best-fit value of all 100 measurements. 

We compa r ed o ur measurements with those of 
iStern fc Laod (|2012[ 1 to cross-check our results with a 
completely independent algorithm (Fig. [2]) and to pro¬ 
vide an estimate of the systematic uncertainties. In con¬ 
trast to our study, their sample contains lower luminos¬ 
ity AGNs at a slightly different redshift range. Thus, 
the samples only overlap partially. For the broad Ha 
FWHM and luminosity, we find relatively tight correla¬ 
tions around the unity relation. The measurements of 
the broad Ha FWHM exhibit a larger scatter and are 



Fig. 2.— Comparison of the fitted FWHM and luminosity of 
the br oad Ha line component between our work and IStern fc La^ 
II2012I1 . for sources in both studies. The dashed red lines show the 
unity correlations. The 1 a scatter about these lines is 0.08 dex 
(left) and 0.04 dex (right), respectively. 



Fig. 3. — Accuracy in recovering the FWHM value of the broad 
Ha line as a function of the signal-to-noise ratio (S/N) of the in¬ 
tegrated flux of the broad Ha line. We use high S/N spectra to 
simulate 250 realizations at various degraded S/N levels. Within 
each S/N bin, we determine the fraction of spectra for which the 
FWHM is within 40% of the original high S/N spectrum (see text 
for details). We show the results for five different objects indicated 
by different symbols. 

slightly skewed toward lower values using our method. 
We inspected those spectra and identify a common char¬ 
acteristic in these cases: the narrow lines are broader and 
weaker than for the other AGNs in the sample. Thus, de¬ 
blending the narrow [Nil]/-Ha from the broad Ha line 
can be challenging. This is particularly true for broad Ha 
lines with a FWHM less than 4000 km s“^. In summary, 
our line-profile fitting method provides robust estimates 
for the FWHM and line flux within their systematic un¬ 
certainties and without the need for human interaction. 

3.2. Estimating Mbh 

The virial black hole mass is estimated by combin¬ 
ing the empirically calibrated BLR size-luminosity re¬ 
lation, derived from reverberation mapping monitoring 
of nearby AGNs, and assuming virialized BLR motions, 
Mbh = /virf^^7?BLRG“^, where v is the FWHM of broad 
lines in the AGN spectrum and /vir is the virial fac¬ 
tor representing the BLR kinematics. Several different 
estimators for Mbh have independently been reported 
in the literature using different BLR size-luminosity re¬ 
lations and virial factors based on different assump¬ 
tions. Here we employ two different calibrations to check 
their potential impact on our results. The first is from 
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iMcLure fc Dunlop (|2004D : 

Mbh / i^sioo y-^V FWHMH^a y 

Mq yW^ergs-^ J V kms-i y 

which is also used bv iShen et alJ (j2009l l to estimate Mbh 
for all unobscured, optically selected SDSS AGN. The 
second calibration is b ased on the most rec ent BLR size- 
luminosity relation of iBentz et alJ (|2009f) . assuming a 
virial fac tor of fyjr = 5.5, w hich is empirically deter¬ 
mined bv lOnken et alJ (|2004ll : 


A^bh 


= 8.13 


(lO^'^ergs 


^5100 


0-52 ^FWHMh/3\' 

kms-1 ) ' ^ ^ 


Throughout the paper, we use the second calibration 
(jBentz et al.l l2009fl for the estimates of the black hole 
masses. In Sect. IQ below we test the robustness of our 
clustering results if the first calibration method is used 
instead. As we demonstrate below, we find very similar 
results using either calibration. 

In order to use our measurements from the Ha line 
to estimate Mbh, we have to replace T 5100 by Thq and 
FWH Mh/ 9 by FWHMHfY in Equations Q and fol- 
lo wing IGreene fc (I2005D (their Equations 1 and 3). 

iDennev et al.l ( 2009ll point out that the inclusion of low 
S/N spectra adds a systematic offset in Mbh estimates 
because the line width is systematically underestimated. 
Therefore, we enforce a lower cutoff in the measured S/N 
level of the broad Ha line flux. We determine the cutoff 
S/N by selecting the five AGN spectra with the high¬ 
est S/N. Then, we artificially degrade the S/N gradually 
down to a continuum S/N of I. For each S/N level, we 
generate 250 spectra and model these individual spectra 
with our line-profile fitting method. 

At an Ha FWHM uncertainty level of greater than 
40%, the measurement uncertainties exceed the com¬ 
monly assumed systematic uncertainties for virial Mbh 
estimates of ~0.3 dex. More than 90% of all of our gener¬ 
ated spectra have FWHM uncertainties of less than 40% 
if we restrict the Ha S/N to values greater than 10 (Fig¬ 
ure El). Gonsequently, we adapt this S/N threshold as a 
limit for selecting objects that yield reliable Mbh esti¬ 
mates. This cutoff leads to a removal of 189 out of 1538 
objects (12% of the sample). We apply this S/N cut and 
compare our measur ed FWHM values to the ones from 
iStern fc Laod (I20I2I1 . The observed scatter in FWHM 
between both methods is 0.08 dex, which corresponds to 
an Mbh error of 0.16. This can be interpreted as a sim¬ 
ple estimate of the minimum systematic uncertainty in 
the Mbh determination, which is reflected in the com¬ 
monly assumed total ~0.3 dex systematic uncertainty of 
this method. 

The main purpose of this paper is to study the clus¬ 
tering of AGN samples with low and high Mbh and 
A/Aedd, respectively. Any uncertainty in the calibration 
of Mbh will affect our entire sample similarly. Thus, our 
results will not depend on the absolute accuracy of the 
Mbh calibration because we are interested only in rela¬ 
tive Mbh values such that we can divide the full sample 
into low and high Mbh subsamples. 


3.3. Estimating L/Ledd 



Fig. 4.— Mbh vs. F/Ledd for luminous, unabsorbed, broad¬ 
line, X-ray selected RASS/SDSS AGNs with a broad Ho luminosity 
of S/N > 10. The color of the symbols encodes the X-ray luminos¬ 
ity of the sources. The lower and right panels show histograms of 
Mbh und T/Tedd, respectively, for the full sample. The dotted 
lines mark the restriction of the parameter space in our analysis to 
remove extreme objects. 

The Eddington ratio A/Aedd is the ratio between 
the bolometric and the Eddington luminosity (Agdd = 
1.26 X IO^®Mbh/M0 erg s“^). The quantity is com¬ 
monly derived by using the inferred Mbh and the optical 
continuum luminosity at 5100 A, adopting a certain bolo¬ 
metric correction factor. Different bolometric correction 
factors used to estimate the bolometr ic luminosity Ahni 
from A.c,inn are quoted in t he literature (iKasoi et al.ll2000t 
iMcLure & Dunloiil 12004 iRichards et al.ll2006aD . While 
these correction factors are subject to uncertainties, they 
are required to estimate the bolometric luminosity if only 
a single-epoch optical spectrum is available. We use a 
correction factor of / = 10.3 (|Richards et al.l [2006all . 
which is consistent with values used by othe r stud ies. 
Again, we use the relation from IGreene fc Hoi (j2005h to 
estimate A 5100 from Lhoc for the computation of Aboi- 

3.4. Defining X-ray Selected AGN Subsamples 

The RASS/SDSS AGNs in our sample do not uni¬ 
formly populate the Mbh - A/Ledd plane (Figure E]). 
Although there is substantial scatter, on average, higher 
A/Aedd are found in AGNs with lower Mbh- While the 
absence of objects in the upper right corner of this plane 
is not an observational bias, but reflects the nonexistence 
of (unabsorbed) RASS/SDSS AGNs with both high Mbh 
and high L/Ledd in the redshift range studied here, the 
lack of objects in the lower left corner reflects an obser¬ 
vational bias. AGNs with low A/Ledd and low Mbh are 
too weak to produce a significant broad Ha line relative 
to the host galaxy continuum. The Ha flux S/N cutoff 
used here primarily removes objects with log {Mbh/Mq) 
~ 7-8. 

Since Mbh and A/Ledd are correlated in our full sam¬ 
ple, and we aim to reveal which of these quantities drives 
the observed X-ray luminosity dependence of AGN clus¬ 
tering, doing a simple cut of the full sample into low and 
high Mbh subsamples, as well as high and low L/Ledd 
subsamples, will not be useful. When creating subsam¬ 
ples that depend on one parameter, the distribution of 
the second parameter in both subsamples must be the 
same. This “matching” of the subsamples is a commonly 
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Fig. 5. — Low (filled down-pointing triangles), medium (filled 
circles), and high (filled up-pointing triangles) Mbh RASS/SDSS 
AGN subsamples with matched L/Ledd distributions. The top 
panel shows the L/Ledd-A/bh plane. The color of the points en¬ 
codes the X-ray luminosity from lowest Lx (blue) to highest Lx 
(red). The solid lines are for visual guidance and mark the di¬ 
vision of the full AGN sample into low, medium, and high Mbh 
subsamples. The bottom and right panels show the L/Ledd ^rid 
Mbh histograms of the low (blue), medium (black), and high (red) 
Mbh subsamples, respectively. Since the medium Mbh subsample 
contains more objects than the low and high subsamples, we nor¬ 
malized the total number of this sample in the L/Ledd histogram 
to match the other samples. For clarity, we slightly shift the x 
values of the low and high Mbh subsamples in the histograms. 

used method in clustering measurements ('e.g.. lCoil et alJ 
I2009H . However, one has to test that the procedure of 
creating matched subsamples is not introducing a bias 
to the clustering results. We will do so below by testing 
differ ent methods to produce matched subsamples (see 
Sect. 15.ip . 

Before creating low and high Mbh subsamples, objects 
that lie outside the range — 2.2 < log(L/LEDD) < — 0.2 
are removed. We then determine the number of objects 
in each L/Ledd bin, using a bin width of 0.2 (loga¬ 
rithmic scale). In each L/Ledd bin, we create three 
subsamples: 30% of objects with the lowest Mbh, 40% 
with medium Mbh, and 30% of objects with the highest 
.^BH (Figure [S]). The procedure creates low, medium, 
and high Mbh AGN subsamples with extremely similar 
(“matched”) L/Ledd distributions but different median 
AIbh (see Table [T]). All Mbh AGN subsamples have me¬ 
dian (log (L/Ledd)) = - 1 . 00 . 

The low (30%), medium (40%), and high (30%) 
L/Ledd RASS/SDSS AGN subsamples with matched 
Mbh distributions are created by first applying a cut 
of 6.7 < log(MBH/M 0 ) < 9.5 to remove extreme sources 
and then following the same approach as described above, 
using bins of Mbh with a bin width of 0.2 (logarithmic 
scale). All L/Ledd AGN subsamples have median (log 
(Mbm/Mq)) = 7.92 - 7.93. 

Figures m and [S] show that divisions into low, medium, 
and high Mbh as well as L/Ledd are very similar to divi¬ 
sions into low, medium, and high Lx- Thus, one cannot 
anticipate if the X-ray luminosity clustering dependence 
is related to L/Ledd or Mbh, as we aim to test here. 

The directly observed parameters FWHM and Lhq do 
not correlate with each other as strongly as Mbh and 
L/Ledd (see Fig. 6 ). However, we decide to use the 
same method to create subsamples defined by FWHM 


and Lho as well. Thus, we again split the full sample 
into three subsamples defined using one parameter, while 
maintaining the same distribution in the other parame¬ 
ter of interest. To create subsamples in FWHM, we first 
limit to objects with 41. 7 < log(LHa/[ergs ^]) < 44.1 
and then use a bin width of 0.2. All FWHM AGN sub¬ 
samples have (log(LHa/[ergs“^])) = 42.72 — 42.73. To 
create subsamples in Lho, we first limit to objects with 
2.9 < log(FWIIM/[kms“^]) < 4.1, and then use bins of 
width 0.1, using the logarithmic values of FWHM. This 
choice of bin width ensures a similar numbers of bins as 
compared to the subsamples defined in Lho • The median 
FWHM of all Lhc subsamples is (FWHM) = 2850-2960 
kms“^. 

As we have to exclude ^13% of the RASS/SDSS AGN 
sample studied in Papers I-HI to fulfill our S/N require¬ 
ment of the broad Ha line flux for reliable Mbh esti¬ 
mates, we create low and high X-ray luminosity subsam¬ 
ples of the reduced AGN sample used here. We use 
a simple 0.1-2.4 keV luminosity cut (with no match¬ 
ing of the distribution of another parameter) of log 
(Lx/[ergs“^]) = 44.29, identical to that used in Papers 
I HI. 

We also apply simple Mi cuts at Mi = —22.10 mag 
and Mi = —22.70 mag to divide the X-ray-selected AGN 
sample into a faint (33% of all objects), medium (34%), 
and luminous (33%) Mi subsamples. Finally, we com¬ 
pute from the SDSS optical spectra the rest-frame ab¬ 
solute magnitude in the band between 5500 and 6800 
A, which includes the Ha line for all objects. We do 
this as the X-ray and optically selected sample covers a 
redshift range of 0.16 < z < 0.36, while the Ha line is 
redshifted outside the SDSS Lband filter above 2 ; ~ 0.28. 
The division into three subsamples of similar size is real¬ 
ized by applying cuts at M|ggQ_gggQ = —21.72 mag and 
-^ 5500-6800 = —22.33 mag. We do not match the distri¬ 
bution of any other parameter for these subsamples. All 
samples used in this paper are presented in Table [T] 

3.5. Defining the Optieally Selected AGN Subsamples 

The optically selected SDSS AGNs (called “quasars” 
in the SD SS literature) ar e draw n from the catalog pro¬ 
vided bv I Schneider et aH (|2010t) Instead of us i ng the 
classic selection in the L-band, iSchneider et al.l (l2010f ) 
use the SDSS f-band, as it is less affected by Galac¬ 
tic absorption. However, this comes with the signif¬ 
icant disadvantage that host-galaxy light might repre- 
sent a significant fract ion of the total flux in the i-filter. 
ISchneider et al.l (|2010[ 1 apply an apparent magnitude cut 
of Mi < —22 and require that objects have at least one 
emission line exceeding an FWHM of 1000 kms“^. Un¬ 
like our RASS/SDSS AGN sample, this sample of 3500 
objects (0.16 < 2 < 0.36) has the footprint of SDSS 
Data Release 7. The DR4-I- and DR7 areas are 5468 deg^ 
and 7670 deg^, respectively, when we consider only the 
area that has a DR7 completeness ratio of /compi >0.8. 
Thus the area occupied by the optically selected AGN 
sample is 1.4 times larger than the area covered by the 
X-ray-selected AGN sample. Additionally, the number 
density (per square degree) of optically selected AGNs is 
1.6 times higher than for the X-ray-selected AGNs. 

We retrieve the individual SDSS spectra of the opti¬ 
cally selected SDSS AGNs and derive the FWHM of the 
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Fig. 6. — Comparison between the X-ray-selected RASS/SDSS AGNs (black lines) and the optically selected SDSS AGN samples (red 
lines). The contours show the location of 10% (dashed lines), 50% (solid lines), and 90% (dotted line) of the full sample, respectively. 
The left figure shows the observed parameters, luminosity and FWHM of the Ho line, while the right figure shows the derived parameters 
L/Ledd Q-nd Mbh- In both figures, a boxcar smoothing with a width of two has been applied. 


broad Ha line and the Lhq through spectral fits (see 
Sect. [S]). This procedure is identical to t he o ne used for 
the RASS/SDSS AGN sample (see Sect. 13.11) . We select 
only objects with S/N > 10 in the flux of the broad 
Ha line. We derive the individual black hole masses and 
L/Ledd for the remaining 2831 AGNs. The X-ray and 
optically selected AGN samples have 807 objects in com¬ 
mon, which is 29% of the total optically selected AGN 
sample. Thus, there is substantial difference between the 
X-ray and optically selected AGN samples. 

Figure[6]compares how the X-ray-selected RASS/SDSS 
AGN sample and the optically selected SDSS AGN sam¬ 
ple span the observed parameter space of Lua and 
FWHM (left panel) and the derived parameter space of 
A/Tedd and Mbh (right panel). In the observed pa¬ 
rameter space, there are two obvious differences between 
these samples. First, the RASS/SDSS AGN sample ex¬ 
tends to lower Ana- This is likely a consequence of re¬ 
jecting AGNs with Mi > —22 in order to exclude ob¬ 
jects with a strong starlight component from the host 
galaxy because an AGN has to have a certain luminos¬ 
ity to outshine its host galaxy. Clearly, X-ray selection 
of AGNs (RASS/SDSS) allows us to extend to fainter 
AGN luminosities, below the optical cut. Thus, at X- 
ray wavelengths we are able to detect AGNs that might 
not outshine their h ost galaxy in the optical (see, e.g., 

I Hopkins et al.l [200^ . At log (AHa/[ergs“^]) > 42.5, 
both samples span a similar parameter space, and the 
dynamic range in FWHM is almost identical in both sam¬ 
ples. 

Differences in the observed parameter space naturally 
translate into differences in the derived Mbh“A/Ledd 
space (Fig. [51 right panel). Compared with the X-ray- 
selected AGN sample, the optical AGN sample has a de¬ 
ficiency of AGNs with low Mbh- The optically selected 
sample also extends marginally to higher Mbh and lower 
A/Aedd- 

We create optically selected AGN subsamples with 
respect to Mbh, A/Ledd, FWHM, and Aho (with 
matched distribution in the other parameter of inter¬ 
est; see Sect. [S3). We have to use slightly different 
limits for some of the parameters because the X-ray 


and optically selected AGN samples cover slightly dif¬ 
ferent parameter spaces. Thus, when we create the 
low, medium, and high Mbh subsamples with matched 
A/Aedd distributions, we use an A/Ledd limit of —2.3 < 
log(A/LEDD) < 0.1. The optical Mbh subsamples have 
median (log (A/Ledd)) = —1-04 to —1.00, while the 
A/Aedd subsamples have (log (Mbm/Mq)) = 8.12 — 
8.13. 

When we create three subsamples in FWHM, we 
limit Ahq to 42.1 < log(AHa/[ergs“^]) < 44.1. The 
FWHM subsamples have {log(AHa/[ergs“^])) = 42.88. 
For the Aho subsamples, we apply a limit of 2.8 < 
log(FWHM/[kms-i]) < 4.2. We find (FWHM) = 
3240 — 3280 kms“^ for the low, medium, and high Lhq 
subsamples. All bin widths are identical to the ones used 
for defining the X-ray-selected AGN subsamples. 

We produce optical AGN subsamples divided by Mi. 
We use simple Mj cuts of Mi = —22.27 and Mi = 
—22.67 to construct faint, medium, and luminous Mi 
subsamples with similar numbers of sources. Finally, we 
create Mglgo-esoo subsamples of the optically selected 
AGN. We use cuts at M||gQ_gggQ = —22.13 mag and 
MItoo- -6800 = —22.53 mag. As with the X-ray-selected 
AGNs, we do not match the distribution of any other 
parameter when creating the subsamples defined by Mi 
and M||gg_gggQ. All sumples are presented in Tabled] 


4. METHODOLOGY 
4.1. Clustering Measurements 

We measur e the two-point correlation function ^(r) 
(lPeeblesll980ll , which measures the excess probability dP 
above a Poisson distribution of finding an object in a vol¬ 
ume element dV at a distance r from another randomly 
chosen object. The AGF measures the spatial clustering 
of objects in the same sample, while the CCF measures 
the clustering of objects in two different samples. We use 
the same approach as described in detail in Section 3 of 
Paper I and Section 4 in Paper HI. Here we explain the 
essential elements of our method. 

We use the correlation estimator of iDavis &: PeeblesI 
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(|1983fl in the form 




DD{r) 

DR{r) 


(3) 


where DD{r) is the number of data-data pairs with a 
separation r, and DR{r) is the number of data-random 
pairs. Both pair counts have been normalized by the 
number density of data and random points. To sep¬ 
arate the effect of redshift distortions, the correlation 
function is measured as a function of two components of 
the separation vector between two objects, that is, one 
perpendicular to (rp) and the other along (tt) the line of 
sight. The parameter ^(rp, tt) is thus extracted by count¬ 
ing pairs on a two-dimensional grid of separations Vp and 
TT. We obtain the projected correlation function Wp(rp) 
by integrating ^(rp,Tr) along the tt direction. 

As in Paper I, we infer the AGN ACF from the CCF be¬ 
tween AGN and the corresponding galaxy tracer set and 
the ACF of the tracer set following ICoil et ^ (1200911 : 


Wp (AGN I AGN) 


[i(;p(AGN|trace)]^ 
Wp (trace I trace) ’ 


(4) 


where Wp(AGN|AGN) and rcp(trace]trace) are the ACFs 
of the AGN and the corresponding tracer set, respec¬ 
tively, and Wp (AGN I trace) is the GGF between the AGN 
and the tracer set (i.e., a galaxy sample). 

The GGF is computed using 

* 41 AGN 4Itrace 

SAGN-trace = 7 :;--t. [O) 

^AGN ^trace 

For our purposes, the use of this simple estimator has 
several major advantages and results in only a marginal 
loss in the S/ N when compared to more advanced esti¬ 
mators (e.g., iLandv fc Szalavl [199311 . The estimator in 
Equation m requires the generation of a random cat¬ 
alog only for the tracer set. Since the random catalog 
should exactly match all observational biases to minimize 
the systematic uncertainties, well-understood selection 
effects are key to creating proper random samples. The 
tracer sets have well-defined selection functions, while 
the AGN samples suffer from selection functions that are 
very complex and difficult to model. 


4.2. Error Analysis 

The measurements in the adjacent bins of the corre¬ 
lation function are not independent. Poisson errors will 
significantly underestimate the uncertainties and should 
not be used for error calculations. Instead, we use the 
jackknife resampling technique to estimate the measure¬ 
ment errors based on the covariance matrix M^-, which 
reflects the degree to which bin i is correlated with bin 
3- 

In our jackknife resampling, we divide the survey area 
into N'T = 100 subsections for the DR4+ geometry (X- 
ray-selected AGN sample), each of which is '^50-60 deg^. 
Since the optically selected AGN sample uses the DR7 
geometry, we divide this area into Vx = 131 subsections 
of roughly equal area. These Vx jackknife-resampled 
correlation functions define the covariance matrix: 


Mij — 


Vx — 1 
Vx 



o.i)- < w{rp^i) > ^ 


X 



< w{rpj) > 


( 6 ) 


We calculate Wp(rp) times, where each jackknife sam¬ 
ple excludes one section and Wk{rp^i) and Wk{rpj) are 
from the kth jackknife samples of the AGN ACF and 
< w{rp,i) >, < w{rpj) > are the averages over all of 
the jackknife samples. The uncertainties represent Icr 
(68.3%) confidence intervals. 

The generation of the covariance matrix for the in¬ 
ferred AGN AGF considers the Nt jackknife-resampled 
correlation functions of the GGF (AGN and correspond¬ 
ing tracer set) and the tracer set AGF. For each jack¬ 
knife sample, we calculate the inferred AGN AGF by 
using Equation m The resulting Nt Wp{rp) jackknife- 
resampled projected correlation functions of the inferred 
AGEs are then used to compute the covariance matrix of 
the inferred AGN AGF. 


5. RESULTS 

We measure high-accuracy GGFs of the different AGN 
samples with the LRG tracer set (see Fig. El). In addition, 
we compute the high-precision AGF of the LRG sample. 
In both cases, we measure rp in the range 0.05-40 h~^ 
Mpc in 15 logarithmic bins, identical to those used in 
Paper III. We compute tt in steps of 5 h~^ Mpc in the 
range tt = 0 — 200 h~^ Mpc. 

To derive 7r„iax, we compute Wp{rp) for a set of TTmax 
ranging from 10 to 160 h~^ Mpc in steps of 10 h~^ Mpc. 
We then fit Wp{rp) over an rp range of 0.3-40 h~^ Mpc 
with a fixed 7 = 1.9 and determine the correlation length 
To for the individual TTmax measurements. As in Papers 
I-III, we find that the LRG AGF saturates at TTmax = 80 
h~^ Mpc. All GGFs saturate at 7r„iax = 40 h~^ Mpc. 
Power-law fits for the AGFs and GGFs are based on 


Wp{rp) = rp 


roY r(l/2)F((7-l)/2) 

r(7/2) 


(7) 


where F(a;) is the Gamma function. 

To derive the clustering properties of the AGN sam¬ 
ples, we follow two different approaches: 


(1) Power-law fits to the inferred AGN ACF: 

We use Eq. E] to infer the AGE for the individual AGN 
sample from the GGF of this sample with the LRG tracer 
set and the AGF of the LRG sample. We fit the data 
points of the inferred AGN AGFs with the expression 
given in Eq. El and derive best-fit ro and 7 values. The 
data are fitted over the range rp = 0.3 — 15 h~^ Mpc 
to be consistent with Papers I-III. Since we measure the 
GGE to infer the AGP, the resulting effective redshift 
distribution for the clustering signal is determined by 
both the redshift distribution of the tracer set and the 
AGN sample: Nccf{z) = Ntrs,cer{z) * Nagn{z). 

In Table [T] we list the redshift range, the median 
effective redshift of Nccf{z) for the corresponding AGN 
samples, the derived best rg and 7 values based on 
power-law hts, and rg for a power-law fit with a fixed 
slope of 7 = 1.9 (for ease of comparison). 

(2) Bias from HOD modeling: In Paper II we develop 
a novel method to infer the HOD of RASS/SDSS AGNs 
directly from the well-constrained GGE of RASS/SDSS 
AGNs with LRGs. In performing the HOD modeling. 
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Fig. 7. — Directly measured cross-correlation functions between AGN subsamples for different parameters with matched distributions in 
the other parameter of interest. In the top row we show the results for the X-ray-selected AGNs, while the bottom row shows results for the 
optically selected AGNs. In each panel, the color encodes the low (red), medium (green), and high (blue) subsample in the corresponding 
parameter. For the sake of clarity, we show only a typical error bar in each panel (black solid vertical line in right upper corner of each 
panel). Note that we show on the y axis Vp x U7(rp) instead of U7(rp) to allow for an easier comparison between the different subsamples 
and to make the values easier to read in the figure. 


we consider that galaxies and AGNs are associated with 
DMHs. A DMH may contain one or more galaxies or 
AGNs that are included in our samples. Using the HOD 
of the LRGs as a template, we constrain the HOD of the 
AGN by fitting the CGF between AGNs and LRGs. A 
few of the GGFs do not have enough pairs on small scales 
for applying statistics. Thus to derive consistent con¬ 
straints for all CGF, we consider only bins with Vp > 0.7 
h~^ Mpc for all CGF fits (as in Paper HI). 

Using the best-fit model of the AGN HODs, we derive 
the bias parameter of the AGN sample by applying 

. _ J bhiMh){NAGN){Mh)4'{Mh)dMh , . 

/(-^AGN)(Afh)<()(Mh)dMh 

where 6h(A7h) is the bias of DMHs with a mass Mh, and 
4>{M'p) is the DMH mass function. In Paper H, we discuss 
different realizations of the HOD modeling depending on 
how central and satellite AGNs in DMHs are included 
or excluded. Here, we apply only “Model A,” assuming 
that, at low redshifts, black hole accretion occurs more 
frequently in galaxies with lower stellar mass than in typ¬ 


ical LRGs. Thus, we assume that (1) the LRG (highest 
stellar mass) occupies the center of the DMH and (2) all 
AGNs are satellites within the same DMH. This is moti¬ 
vated by observations of AGN “downsiz ing” (lUeda et al.l 
I2003I[Hasingedl^OOSHMivaii et al.ll2015h . 

As we are interested here (and in Paper HI) only in 
the large-scale bias values and the estimate of the typical 
halo mass, our results do not depend significantly on the 
choice of the one-halo term model (model A, B, or C). 
For example, in the case of the total RASS AGN sample, 
the best-fit HOD-derived bias values change by less than 
0.02 among models A, B, and C. This is a minor fraction 
of the statistical uncertainty in the bias measurement 
(±0.09). 

The bias parameters derived by power-law fits and 
HOD modeling can deviate significantly. As discussed in 
Paper HI (see Section 5.4. of Paper HI for more details), 
HOD-derived bias parameters are strongly preferred over 
those from power-law fits for various reasons. 

Inconsistencies between the HOD-derived b{z) and the 
results from power-law fits are caused by (1) obtaining 
power-law fits to the inferred AGN ACFs, which have 
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TABLE 1 

Properties of Inferred AGN ACFs and Derived Quantities 


AGN Sample 

Number 

Median 

10th,50th,90th 

ro 

7 

To,7=1.9 

b{z) 

log 

Name 

of Objects 

^efF 

Percentile 

(h-i Mpc) 


(h-i Mpc) 

HOD 

(h-l Mq) 


X-ray 

Selected AGN — RASS/SDSS 

AGN — SDSS Data Release 4+ 



total RASS AGN 

1349 

0.27 

43.70,44.17,44.68 

4 02'’"’^'^^ 

^•^^-0.57 

1 

l-O^_0.14 

^.u»_o 39 

1 q9 + 0.11 

i.oz_0 09 

13.20l°;« 

low Lx RASS 

858 

0.25 

43.63,43.99,44.23 

o 1 9+0.57 

1 f^p;+0.25 

i.OO_Q 21 

q 9Q + 0.41 
o.^o_o 46 

1 99 + 0.14 
-^•^^-0.12 

13.08l°l® 

high Lx RASS 

491 

0.29 

44.34,44.53,44.87 

r 00 + 0.68 

0.00_Q gy 

J--oo_ 0 _i 8 

5-4+o:?3 

1 5q+0.i5 
i.09_Q 17 

13.4+0;« 

low Mbm RASS 

410 

0.22 

7.07,7.54,8.01 

9 -T-i +0.53 

2 07+1-46 

2.65t°;56 


12 50+0-43 

1 Z.DU_q 39 

medium Mbm RASS 

525 

0.27 

7.55,7.92,8.36 

4 Q4+0.64 

1 -Ol_ 0 .i 8 

4 39+0.52 
4.oz_0.58 

1 OO+0.15 
l.dO_Q 17 

1 97+0.16 
10.Z1_Q 23 

high Mbm RASS 

400 

0.29 

7.92,8.38,8.91 

4-60t'i’;n 

i.OO_Q 22 

4 57+0.78 

4.D(_q 92 

1 -O^_ 0 .l 6 

13.49t°;i« 

low -L/Ledd RASS 

418 

0.23 

-1.79,-1.36,-0.98 


1 45+0.26 
l-^0_o.21 

q 1 7+0.54 
'^•1' -0.64 

1 30+0.18 

l-oU-0.18 

13.2i10;21 

medium -L/Ledd RASS 

522 

0.27 

-1.36,-1.00,-0.71 

2.87t\fo 

1 50+0.40 
i.ou_o.33 

2.89t°;« 

1 29"*”01^ 
i-^y-0.15 

1315 +0.16 
19.1D_q 23 

high L/Ledd RASS 

404 

0.29 

-1.01,-0.68,-0.35 

4.c>o_q 79 

2 14+0-22 
^•1^-0. 19 

3-92lg:^^ 

1 18+0.15 
1-1O-0.14 

12.97t°;^2 

low FWHM RASS 

415 

0.27 

1280,1910,2470 

2-69l°jJ 

Z.OU_Q 4Q 

9 q9 + 0.71 

-i-o^-1 01 

1 n9~^o.i2 
i.uy_o.i2 

12.84l°l° 

medium FWHM RASS 

525 

0.26 

2340,2890,3690 

2-9ot^:i 

1 ok+0.24 
i.00_Q 28 

3 7Q+0.64 

0. /U_o.76 

1 47+0.17 

l-^'-0.18 

13.38l°is 

high FWHM RASS 

404 

0.27 

3650,4840,7920 

o K9+0.92 

1 ro + 0.25 
i.00_Q 26 

q 09 + 0.65 

o.oz_Q _77 

1 53+0.17 

i.OO_Q ig 

1 3 43+0.15 
19.4C>_q 20 

low Lboc RASS 

414 

0.22 

41.99,42.29,42.53 

2 14+0'96 

1 39+0.56 
i.oy_o 32 

2 5Q+*^-^4 

z.uy_o 66 

1 91 +0.18 

-^•^-^-0.19 

13 10+0-23 

19.1U_q 35 

medium L^a RASS 

526 

0.27 

42.52,42.74,42.96 

J--yo_1.97 

1 QQ + l^^ 

i.oo_o 36 

9 -I q+0.64 

71 

1 *^1 +0.16 
^•'^^-O.IO 

13 ic+0.19 
19.1©_Q 24 

high Lila RASS 

403 

0.29 

42.93,43.22,43.63 

o K9+O.92 
'^•^^-1.57 

1 cro + 0.25 

1 .oo_Q 26 

q 09 + 0.65 

o.oz_o_78 

1 35+0.18 

1-'JO_0.16 

1 90+0.19 

±o.^o_Q 23 

faint Mi RASS 

444 

0.23 

-21.26,-21.80,-22.04 

3.21«;69 

1 59+0.69 
i.oy_Q 42 

q q9+0.58 

9.9Z_o 70 

1 10"^0.26 
1-11^-0.14 

12.901°;39 

medium Mi RASS 

460 

0.27 

-22.15,-22.39,-22.63 


1 71 +0.36 
' -^-0.29 

3 55+0.61 

9.00_q 72 

1 91 +0.i9 

i-^i-0.16 

13.04tO;25 

luminous Mi RASS 

445 

0.29 

-22.76,-23.13,-23.96 


-*-•^'-0.18 

2 

-‘•i^^_0.65 

1-40-0.16 

13 35+0-16 

19.90_q 19 

faint Mlfptj.ggpo RASS 

444 

0.22 

-20.84,-21.39,-21.66 

SAStiil 

4- /y_o.4i 

3 1 9+0.59 
o.iy_o.7i 

1 1 q+0.25 
1-1O-0.19 

12 97+0.35 

iz.y (_0 42 

med. MIfoVesoo RASS 

460 

0.27 

-21.79,-22.02,-22.27 

2 44"*"^-^^ 

1 35+0.29 

i.oo_Q 28 

3 09+0.57 
o.ou_Q 67 

1 95+0.16 

l-^0_Q 17 

13 1 o+O-^O 
19.1 U_q 28 

lum. Mlfptj.ggoo RASS 

445 

0.30 

-22.42,-22.81,-23.67 

00 

1 70+0-15 
l + U_o.i 6 

4.^U_q 79 

1 44+0.16 

l-^^-0.17 

13 31 +0.15 
19.9i_Q 21 


Optically Selected AGN — SDSS AGN — SDSS Data Release 7 



total SDSS AGN 

2831 

0.29 

-22.08,-22.45,-23.34 

-0.28 

T 80 IH 9 

5 06~*"*^'^^ 
O.UD_q 26 

1 q9 + 0.07 

i.oz_9 Qg 

1 q 1 q+0.08 
19.18_q Q9 

faint Mi SDSS 

932 

0.28 

-22.03,-22.13,-22.24 

5 S4+0-42 
^•'^^-0.51 

i.801°;2° 

5 34+0.40 

0-04_9 43 

1 47+0.09 

1-4^'-0.14 

13.361^6 

medium Mi SDSS 

966 

0.29 

-22.30,-22.45,-22.63 

3 70+0.44 
'J- ' ^-0.59 

1 70+0.16 

4-/0_o 16 

q Q7+0.36 
' -0.39 

1 20 ’^^’^^ 
i-^y_o.io 

13.00l°;i^ 

luminous Mi SDSS 

933 

0.30 

-22.73,-23.06,-23.86 

^•^0-o_64 

1 53+0.12 

i.OO_o 12 

^•0^-0.45 

J^-4U_o.i2 

13 26+^+^ 
19.ZO_q 15 

low Mbm SDSS 

857 

0.27 

7.37,7.82,8.35 

4.^4_o 56 

1.7+0;15 

4-26tn5 

1 90+0.12 

1-^0-013 

13.14l0;i5 

medium Mbm SDSS 

1114 

0.30 

7.70,8.08,8.59 

4 07 + 0.40 
4.y/_o.48 

l-98t°;i^ 


1 1 o+0.13 
1-1O-0.09 

12.96l°;i® 

high Mbm SDSS 

845 

0.31 

8.04,8.45,8.97 

A 79+0.63 

4. ^^_o.82 

i.OU_o.i 2 

4 75+0.43 

^+O_o.47 

1 fiO'^o.io 

1-OU_0.08 

13 45+0.09 
19-40_q Q 7 

low L/Ledd SDSS 

859 

0.28 

-1.84,-1.31,-0.88 

3 40+0-55 

i t::o + 0.14 
^•O»_0.15 

'5-0^-0.44 

1 21 '^ 0 .ll 
i-^i-o.io 

13.03l°;J® 

medium LILbbu SDSS 

1113 

0.30 

-1.53,-1.04,-0.67 


i-o^_0.l3 

'^-oU-o 46 

1 40"^0.11 
1-4U-0 11 

13.26l°i2 

high L/Ledd SDSS 

845 

0.31 

-1.17,-0.74,-0.40 


i-o^-0.15 

A 7q+0.50 
^•''^-0.55 

1 49+0.15 
J^-4y_o.i2 

13.35l0;i4 

low FWHM SDSS 

857 

0.30 

1490,2110,2560 

4 18+0-49 
^•-*^^-0.56 

1 92"^'^4L® 

i.yz_0 19 

4 17+0.49 
^•l'-0.55 

1 11+0.13 
1-11-0.08 

12.84l®li 

medium FWHM SDSS 

1120 

0.29 

2680,3260,4040 

K qo+0.43 

1 75+0.i3 
l-'0_o.i3 

5 1 9+0.38 
O-l^-0.40 

1 41+0.10 

1-41-0.09 

13.281°+ 

high FWHM SDSS 

847 

0.29 

4300,5500,8460 

63 

1 OR+0-14 

l-bo_0,i4 

4 71+0.42 

^•'1-0.46 

1 44+0-10 

i-4^^-0.12 

1 SI q9 + 0 .i 0 
±o.oz_Q 15 

low Lila SDSS 

859 

0.27 

42.40,42.61,42.71 

0 79 + 0.51 

0. <^_o.72 

1 53+0.16 
i.DO_o.i7 

3 79+0.41 
'5- /y-0.45 

1 95+0.13 

13.10l°+ 

medium Luoc SDSS 

1119 

0.30 

42.76,42.89,43.04 

^•^'J_ 0.68 

1 59"^^1^ 

i-oy_o.i 6 

4 47+0.39 
4-41-0 42 

1 q9 + 0.12 

i.OZ_Q IQ 

13.16t+t 

high Lila SDSS 

846 

0.31 

43.07,43.26,43.63 

^•o^_0.80 

l-00_o.i4 

4.0U_o.5i 

1 59+0.i 3 

i.oz_o 13 

13 . 381+2 

faint M5|gg_gggQ SDSS 

933 

0.27 

-21.70,-21.97,-22.10 

5.96l°;« 

1 ss+0.16 
l-oo_0.l4 

5 95+0.40 

o.yo_Q 43 

1 50+0.09 

i.OU_o 14 

13.401°;°° 

med. MlfgVesoo SDSS 

964 

0.30 

-22.17,-22.31,-22.49 


1 55+0.20 
i.OO_Q 22 

3 95+0.40 

o.yD_Q 45 

1 - 14 + 0.10 
l-l^-O.ll 

12.891+® 

lum. M|fg‘g_gggg SDSS 

934 

0.31 

-22.59,-22.89,-23.66 


1 594"0.13 

i-oy_o.i3 

4 55+0-48 

4.0D_q 53 

1 45+0-12 
1-40-0.13 

1 q9 + 0.12 

19 . 9 Z_q 12 


Note. — Samples displayed in italics are those that are created to match distributions in the other parameter of interest 
(e.g, samples split in Mbm have matched distributions in L/Ledd and vice versa; samples split in FWHM have matched 
distributions in Lna). The fourth column lists either the log (Lx/[ergs“^]), log {Mbm / Mq) , log (L/Ledd), FWHM (kms“^), 
log (LHa/[ergs“^]), absolute magnitude in the i-band (Mi in mag), or absolute rest-frame magnitude in the 5500-6800 A 
band (M 5500 _ 68 oo in iTia-g) of the studied sample. In this column, the 50th percentile corresponds to the median value. For 
the total sample, we state in this column the Lx or the Mi value. All samples cover a redshift range of 0.16 < 2 ; < 0.36 and 
contain only objects with S/N > 10 in Ha flux. Values of ro, 7 , and ro.T-^i.g are obtained from a power-law fit to Wp(Tp) over 
the range rp = 0.3-15 h~^ Mpc, using the full error covariance matrix and minimizing the correlated values. 
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much larger uncertainties than the directly measured 
CCFs between AGNs and LRGs (the HOD uses GGFs 
instead of the AGN AGFs), and (2) fitting the inferred 
AGN ACF in a range Vp = 0.3-15 h~^ Mpc. The latter is 
also applied by other studies because the clustering signal 
is not well constrained above Vp ~ 15 h~^ Mpc. How¬ 
ever, such a fit includes scales in the nonlinear region 
in which the bias-DMH mass relation, based on linear 
theory, should not be applied. We decide to do so never¬ 
theless to be consistent with Papers 1 and HI, as well to 
allow the reader to perform a direct comparison to other 
studies. 

Based on the HOD-derived bias parameter {h{z) 
HOD), we derive the typical DMH mass occupied by 
the A GN sample. Using Equation ( 8 ) of iSheth et al.l 
(1200111 and the im proved fit for this equation given by 
iTinker et al.](l200^ . we compute the expected large-scale 
Eulerian bias factor for different DMH masses at differ¬ 
ent redshifts. Comparing the observed b value from HOD 
modeling with the DMH bias factor from ACDM cosmo¬ 
logical simulations provides an estimate of the typical 
DMH mass (&DMH(AfQ^j^jj) = &obs.hod) in which the 
different AGN samples reside, as listed in the last col¬ 
umn of Table [TJ Small differences in the median effective 
redshift between different subsamples do not significantly 
change the derived typical DMH mass. The low and high 
Mbh samples of the X-ray-selected AGN sample have 
the largest difference in median effective redshift when 
we consider only samples that have matched distribu¬ 
tions in the other parameter of interest. Computing the 
typical DMH mass at these different redshifts leads to a 
difference of only Alog(MQ^jj) = 0.08 (e.g., for a bias 
value of 1.20, log 22 / 1 ^”^ ^q])= 12.83 and 

log (AfB^jj,z=o. 3 o/[^ ^ Af 0 ])= 12.75). 

We emphasize that should not be compared 

between different studies. Various conversions are used in 
the literature to derive from 6 dmh- In particular, 

for the same bias parameter, studies based on optically 


selected AGN samples derive Af^^jj values up to ^0.6 
dex lower than studies using X-ray-selected AGN samples 
(depending on the actual bias factor and r edshift of the 
sampl e). The use of the improved fi t from I Tinker et al.l 
(j2005D compared to Equation ( 8 ) of ISheth et al.l ( ~001h 


alone can account for an Af^^jj difference of ^0.3 dex. 
Thus, instead of blindly using the derived values 

from different studies, one should recompute these values 
from the bias parameters in the same manner. 

For the X-ray and optically selected AGN samples, we 
find a clustering dependence on AIbh, in that subsamples 
with high Mbh cluster more strongly than their low Mbh 
counterparts. No significant clustering dependence on 
L/Ledd is observed. We also detect a weak FWHMna 
clustering dependence, in that AGN with low FWHMna 
are less clustered than their high FWHMho counterparts. 
For the X-ray AGN sample, which has a larger dynamical 
range in L^a than the optical sample, we do not find a 
significant dependence on Lho- 


5.1. Robustness of the Clustering Measurements 

In this section we verify the reliability of our cluster¬ 
ing results by altering the selection of the different AGN 
subsamples with respect to different parameters. Mod¬ 


erate changes should not influence the clustering results 
significantly. Therefore this should serve as a test of the 
robustness of our results. 

In Section 13.21 we clarify why we do not use black 
hole mass estimates from objects with Ha luminosity 
S/N < 10. For the first check, we lower this thresh¬ 
old to S/N=5, generate the various AGN subsamples, 
and compute their clustering. Comparing the subsam¬ 
ples with respect to the same parameter (e.g.. X-ray lu¬ 
minosity, Mbh, T/Ledd), the newly derived constraints 
agree well within their la uncertainties with the findings 
of the S/N=10 samples. The same applies when we apply 
a very conservative threshold of S/N=20. 

The intrinsic scatter of approximately 0.3 dex in the 
SMBH mass estimate could have a substantial effect on 
the observed correlations between Mdmh and Mbh as 
well as L/Ledd- In Section 1^^ we use simulations to 
explore the impact of SMBH mass errors on the cluster¬ 
ing correlations and show that the expected modification 
of the correlations is very mild. 

In this section, we investigate possible effects of the 
scatter in the SMBH mass estimates from the observa¬ 
tional point of view. Two different approaches for esti¬ 
mating M bh b ased on the optical spectra were discussed 
in Section !??^ Throughout the paper, we use the rela¬ 
tion given in iBentz et al.l (I2009H fsee Eg. Her e we 
recompute Muu based on lMcLure fc Dunlool (|2004[1 (see 
Eq. [T]). For the X-ray and optically selected AGN sam¬ 
ples, t he Mbh values estimated from lMcLnre fc Dnnlod 
(l2004li are on average 0.2 dex smaller than the ones 
from IBentz et al.l (120091) . We repeat the clustering mea¬ 
surement of low and high Mbh and L/Ledd AGN sub¬ 
samples, respectively. We find trends identical to those 
reported in Table H] The result s for all subsamples 
based on lMcLure fc Dunlod (120041 ) agree well wi t h thei r 
corresponding subsamples based on IBentz et al.l (j2009t) . 
within the la uncertainties. Thus, different estimators 
for the SMBH mass, which are subject to different sys¬ 
tematic errors, give extremely similar results. We con¬ 
clude that, in our study, uncertainties related to the clus¬ 
tering measurement itself (i.e., the moderate number of 
AGNs) dominate over systematic SMBH mass errors. 

A key ingredient of our paper is that we split the AGN 
sample in such a way that we match the distribution of 
another parameter of interest (see Sect. 13.41) . To test this, 
we alter the procedure of creating the different subsam¬ 
ples. First, we fit a regression line to the data and split 
the sample into objects above and below the regression 
line. The resulting subsamples have moderately deviat¬ 
ing distributions in the other parameter of interest. This 
approach yields consistent clustering results with respect 
to our final results presented in Table [T] 

Second, we repeat the generation of low, medium, and 
high Mbh and L/Ledd samples with matched distribu¬ 
tions in the other parameter of interest. This time, we 
split the sample into 40% (low), 20% (medium), and 40% 
(high). All corresponding samples agree with the original 
samples, within the combined la uncertainties. Except 
for the optical L/Ledd subsamples, we also find identical 
trends for all other parameters in the X-ray and optically 
selected AGN samples. For the optical L/Ledd subsam¬ 
ples, the split into 30% (low), 40% (medium), and 30% 
(high) suggests a weak and tentative positive L/Ledd 
clustering dependence (^l.Tcr between the low and high 
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L/L-e^Y) subsamples; see Tabled]). When we split the 
same sample into 40% (low), 20% (medium), and 40% 
(high), we find &iow = 1.351°([J, ^medium = 1.33l°lg, 
and 6 high = I-SSIqIi- Thus, we find no evidence for an 
A/Tedd clustering dependence. 

Third, we generate subsamples that have perfectly 
identical distributions in the other parameter of interest 
(instead of very similar distributions). To do so, we com¬ 
pute in each bin which of the subsamples (low, medium, 
or high) contains the most objects. We then randomly 
select and add again objects from the subsamples with 
the lower number of objects until all three subsamples 
contain an equal number of objects in this bin. In other 
words, we scale up the number of objects in a subsam¬ 
ple to that of the subsample with the highest number of 
objects. This method allows multiple entries of single ob¬ 
jects in a given subsample. Ultimately, in each bin only 
one or two objects are added to subsamples. Thus, these 
subsamples are extremely similar to the matched sam¬ 
ples that we use throughout the paper. Consequently, 
it is not surprising that both methods result in virtually 
identical clustering signals. 

In Fig. [5| (right panel histogram) the low, medium, 
and high Mbh subsamples overlap considerably. As 
a fourth test, we now create low and high Mbh and 
A/Aedd subsamples that have completely separated dis¬ 
tributions, but still match in the L/Ledd and Mbh dis¬ 
tribution, respectively. We divide the X-ray and optical 
AGN samples at Mbh = 10® Mq and log (A/Ledd) = 
— 1. We weight the samples by adding sources in a 
subsample multiple times until a matched distribution 
is created. One major disadvantage of this approach 
is that we have to limit substantially the parameter 
space of the other parameter to match its distribution 
in both subsamples (see Figs. |31[5]). For the X-ray- 
selected RASS/SDSS AGN samples, we find 6 iowMbh = 

1 rkr+0.18 t _ -| r)^ + 0.21 7 -« -|q + 0.18 

i.UO_Q 2 ^0, OhighMBH ~ -^-^^-0.20’ ^IowL/Ledd ~ -‘-■-*-'^-0.18’ 

and 5highL/LEDD = 1-17^0 i 4 - The bias values for the 
optically selected SDSS AGN samples are 6 iowMbh = 

1 n(^+0.14 7 _ 1 OO+0.16 L _ 1 on+O.lO 

OhighMBH “ ^■'40-0.15’ "IowL/Ledd ~ ^■'^^-0.137 

and fehighL/LEDD = l- 28 io' 2 i- Due to significant restric¬ 
tion of the parameter space and thus a substantial de¬ 
crease in the number of AGN, this approach leads to less 
significant results. 

Fifth, we combine the X-ray and optically selected 
AGN samples. We divide the combined sample into the 
highest 30%, middle 40%, and lowest 30% with respect to 
the corresponding parameter. We also combine the mid¬ 
dle and lowest samples (70%). Again, the different sub¬ 
samples have matched distributions in the other parame¬ 
ter of interest. We find a 2.6(T difference when comparing 
the AGNs with the 30% highest to the 70% lowest SMBH 

masses (&70%1 owMbh “ 1T9_q gg, &highMBH “ 2o)- 

Last, we drop the requirement of having subsamples 
with matched distributions in the other parameter of in¬ 
terest. We naively create samples of low, medium, and 
high Mbh, A/Ledd, FWHM, and Ahq simply by divid¬ 
ing the samples into the highest 30%, medium 40%, and 
lowest 30% with respect to the corresponding parameter. 
Dividing the sample according to Mbh, for example, in 
this manner leads to subsamples with very different me¬ 
dian A/Ledd values. We present the results for subsam- 


TABLE 2 

Clustering Properties with Respect to A^bh and A/Ledd 
When Using Unmatched Distributions in the Other 
Parameter of Interest 


AGN Sample 

Median log 

Median log 

b(z) 

Name 

(AtBH/At©) 

(A/AedD j 

HOD 

X-ray-selected RASS/SDSS AGN 

- Data Release 4+ 

low Mbm RASS 

7.49 

-0.76 

1 

-L'UU_o.i5 

medium Mbm R-ASS 

7.92 

-1.01 

-L-oO_o 18 

high Mbm RASS 

8.46 

-1.33 

1 40+0-15 

-L-^y_o.i6 

low L/Ledd RASS 

8.29 

-1.42 

1 41+0.17 

l'^l-0.16 

medium L/Ledd RASS 

7.92 

-1.00 

1 4Q+0.20 

high L/Lbdd RASS 

7.62 

-0.64 

i.uy o.i2 

Optically Selected AGN 

- SDSS AGN 

— SDSS Data Release 7 

low Mbm SDSS 

7.72 

-0.72 

1 OQ-pU.ii 

l-^o_0.l6 

medium Mbm SDSS 

8.12 

-1.04 

1 ik+0.08 

J^-J-©_0.07 

high Mbm SDSS 

8.58 

-1.42 

1 e:Q+0.07 

T.U0_Q QQ 

low L/Abdd SDSS 

8.51 

-1.46 

1 43+0.11 

medium A/Abdd SDSS 

8.08 

-1.03 

1 42+^+^ 

high A/Lbdd SDSS 

7.76 

-0.65 

1 19+0.15 


Note. — An explanation of the columns is given in 
Table [TJ 


pies in Mbh and L/Ledd with unmatched distributions 
in the other parameter (for X-ray and optical AGN) in 
Table [2| Despite the “mixing” of the two parameters, 
Mbh and L/Ledd, the trends seen in the HOD-based 
bias parameters are similar, though not identical, to the 
subsamples with matched distributions. We will discuss 
these results below in Sect. ED in more detail. 

As described above, we consider only data bins for all 
CCFs with Pp > 0.7 h~^ Mpc for the computation of the 
HOD-based bias parameter. As a last test, we recom¬ 
pute the HOD bias parameter using bins with Pp > 0.5 
h~^ Mpc. Again, the results agree well within the Icr 
uncertainties. The few exceptions are caused by a very 
low data point at Vp ^ 0.6 h~^ Mpc. In all of these cases, 
there are at most 10 pairs that contribute to this data 
point. Thus, statistics, as mentioned before, should 
not be applied. The combination of the various tests 
listed here provides convincing evidence that our results 
are not significantly influenced by systematic effects and 
demonstrates their robustness with respect to moderate 
changes in our methodology and sample selection. 

5.2. Exploring the One-halo Term Dependence of 
Different A GN Parameters 

The HOD modeling introduced in Paper H of this se¬ 
ries allows us to explore properties of the one-halo term, 
that is, at small separations where two objects occupy 
the same DMH. More importantly, we are able to provide 
the full distribution of the number of AGNs as a func¬ 
tion of Mdmh, instead of quoting only a typical value. 
As a parameterization of the AGN HOD, we use a sim¬ 
ple truncated power law, assuming that all AGNs are 
in satellites (’Model A’ in Paper H, for more details see 
Paper H): 

(A^agn.s) oc Mdmh“ ©(Mdmh — Mcr) (9) 

where Q{x) is the step function (=1 at a: > 0; =0 at 
a; < 0), Mcr is a critical (minimum) DMH mass below 
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log M„ 

Fig. 8.— Confidence contours of the HOD-derived parameters a and log Mcr for the low and high Lx, Mbh, -^^/^edd, FWHM, and 
RASS/SDSS (X-ray-selected) luminous, broad-line AGN subsamples (from left to right). The blue contours show the high AGN 
subsamples in one parameter, and the red contours show the low subsamples. The confidence intervals correspond to the = 1-0 and 
2.3 levels. The green dashed lines show the mean DMH mass in log (M^mh/^”^ ^©1) derived from the model parameters. Every other 
contour level is labeled. 



log 

Fig. 9. — Similar to Fig.|8] here showing the confidence contours of the HOD parameters a and log Mcr for the low and high Mi, Mbhi 
^/^BDD) FWHM, and Lhq optically selected SDSS AGN samples (from left to right). 


which the HOD is zero (where the DMH does not contain 
an AGN), a is the power-law slope of the HOD above 
Mcr, and Vagn,s the number of (satellite) AGNs in the 
same DMH. 

While a model with an adequate mix of central and 
satellite AGNs is more realistic (models B and C in Pa¬ 
per H), we here show the results of our simplest model 
because the purpose of this section is to highlight the dif¬ 
ferences between samples, including the constraints from 
both the one- and two-halo terms, while the constraints 
on bias values come only from the two-halo term. 

To achieve adequate constraints, we divide the full 
AGN sample into only two subsamples (instead of three, 
as above) for each parameter of interest. We follow the 
description of Sect. 13.41 and create low and high sam¬ 
ples for Mbh, T/Tedd, FWHM, and Lua with matched 


distributions in L/Ledd, A^bh, Lua, and FWHM, re¬ 
spectively. We list their clustering properties in Ta¬ 
ble El The low and high Lx and Mi samples use a 
simple cut (log (Lo.i-2.4kev/[ergs“© = 44.29, Mi = 
—22.4). For the faint Mi sample {{Mi) = —22.18), we 
find fefaint = 1.45to ) 2 , while the luminous Mi sample 
{{Mi) = —22.79) yields hum. = l-29to'io (l-2cr differ¬ 
ence). The Mi clust ering dependence will be discussed 
in detail in Sect. 16.31 

We run the HOD modeling for all low and high AGN 
subsamples using bins with rp > 0.5 h~^ Mpc. To do 
so, we assume that statistics can be applied if a data 
point contains at least 15 pairs. Lowering the rp limit al¬ 
lows us to explore the one-halo term via HOD modeling 
in more detail. We derive the best fit for all of the low 
and high parameters individually and show their results 
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TABLE 3 

Clustering Properties with Respect to Mbh and L/Ledd for 
Divisions into Two Subsamples with Matched Distributions 


AGN Sample 

Median 

Median log 

Median log 

b(z) 

Name 

^efF 

(Mbh/Mq) 

(L/Ledd) 

HOD 

X-ray selected AGN - RASS/SDSS AGN - SDSS Data Release 4+ 

low AfBM 

0.24 

7.66 

- 1.00 

1 -|4+u.ic} 

high Mbm 

0.29 

8.23 

- 1.00 

1 42+0.15 

low L/Ledd 

0.25 

7.93 

-1.24 

1 oq+0.15 

l-oU_o.i3 

high L/Ledd 

0.28 

7.92 

-0.77 

1 9Q”t"0.13 



Median 

Median log 




FWHM 

Lho 


low FWHM 

0.26 

2170 

42.72 

1 1 c+u.iy 

high FWHM 

0.27 

4050 

42.72 

1 47+0.13 

low Lhq 

0.24 

2890 

42.45 

1 -,^.+0.16 

1-I0_o.i4 

high Lua 

0.29 

2890 

43.04 

1 44+0.13 

Optically Selected AGIN 

— SUSS AGJN 

- SUSS Uata Helease 7 

low Mbm 

0.28 

7.89 

- 1.02 

1 oQ+o-oy 

high Mbm 

0.30 

8.32 

-1.04 

1 4 Q+ 0.10 

low LjLBou 

0.28 

8.12 

- 1.22 

1 qO + 0.08 
i.oo_Q 09 

high L/Ledd 

0.30 

8.12 

-0.83 

-1 00+0.09 

l.o»o.ii 



Median 

Median log 




FWHM 

LHQ! 


low FWHM 

0.26 

2440 

42.88 

1 97+u.iU 
- 0.08 

high FWHM 

0.27 

4550 

42.88 

1 45+0.08 

i-^o_0 11 

low Lhc 

0.24 

3250 

42.69 

1 oq+0.09 
l-oU_0 j0 

high Lua 

0.29 

3280 

43.10 

1 42+0 09 
l-+^0.09 


Note. — An explanation of the columns is given in Ta¬ 
ble [TJ The units of FWHM are kms“^ and log Lna are log 
(.bHa/[ergs“^]). Note that the units and column descrip¬ 
tions of the optically selected SDSS AGNs are identical to 
the ones used for the X-ray-selected RASS/SDSS AGNs. 

in Figs. | 8 ] and [9] These two-dimensional plots show a 
clearer difference between the subsamples than collaps¬ 
ing the information into one dimension such as the bias 
parameter. 

All low and high subsamples (independent of the stud¬ 
ied parameter) are consistent with a < 1. As pointed 
out in Paper II, HOD analyses of galaxies over a wide 
range of absolute magnitudes and redshifts find a ~ 
1.0—1.2. Gomparing the AGN and galaxy results implies 
that models are preferred in which the fraction of satel¬ 
lite AGNs to satellite galaxies decreases with increasing 
DMH mass. Fig.[5]shows evidence that high FWHM and 
broad-line Lna AGNs show an even steeper decrease with 
DMH mass (i.e., lower a). 

The confidence contours of the optically selected 
AGN subsamples (Fig. [9|) are narrower than for the 
X-ray-selected sample. This is due to the larger number 
of optical AGNs than X-ray-selected AGNs. However, 
the general two-dimensional appearance of the confi¬ 
dence contours are very similar between the X-ray and 
optically selected AGN samples. We note that, for the 
Mbh subsamples, the low and high optically selected 
AGNs have different a and Mcr values, while the low 
and high Mbh X-ray-selected AGN subsamples have 
consistent a and different Mcr values. 


6 . DISCUSSION 

In this section, we use our results to discuss the physi¬ 
cal origin of the observed Lx clustering dependence. We 
also compare our results to other studies and discuss why 
we do not detect an Mi clustering dependence when the 
clustering results between the X-ray and optical AGN 
samples are so similar. Finally, we compare our results 
to predictions from state-of-the-art semianalytic cosmo¬ 
logical simulations. 

6 . 1 . Mbh o,s the Origin of the Lx Clustering 
Dependence 

We first verify that the weak X-ray luminosity depen¬ 
dence of the clustering strength found in Paper I is still 
present in the reduced AGN sample when we exclude 
^13% of all objects to reliably estimate Mbh- The de¬ 
pendence is detected at the 1.4cr level when considering 
“ 6 ( 2 ;) HOD” in Table [T] (and is 2.5 ct comparing the rp 
values at fixed 7 ). Using the full X-ray sample, we re¬ 
ported a l.Scr detection in Paper H using HOD-derived 
bias values. 

We split the X-ray-selected AGN sample into three 
subsamples with respect to Mbh (with matched L/Ledd 
distributions) and show the results in Fig. (TU] (left) as a 
blue line. The X-ray sample shows a steady increase of 
the bias parameter with Mbh- The difference between 
the lowest and highest Mbh subsamples in the X-ray- 
selected AGN sample is 2.3tT. 

For the X-ray and optically selected AGN sample (red 
line), we detect the highest clustering in the subsample 
that contains the 30% most massive black holes. For 
both AGN populations, we combine the low and the 
medium subsamples and compute their combined clus¬ 
tering strength. As a result of the increased sample size, 
the uncertainties decrease. The bias difference compared 
to the high Mbh subsample is 2 .Oct and 2.7 ct for the X-ray 
and optical AGN samples, respectively. We show these 
data points at the median Mbh values of each subsample 
in Fig. [ini 

The X-ray AGN sample shows no hints of a correlation 
between clustering strength and L/Ledd (Fig. [TUI right 
panel). The optical AGN sample has a slight positive cor¬ 
relation in that AGNs with higher L/Ledd cluster are 
slightly more clustered than their lower L/ Ledd counter¬ 
parts. However, this ~1.7 ct difference disappears when 
we split the sample into 40% (low), 20% (medium), and 
40% (high) bins in L/Ledd (see Fig. [TUI right panel). 
In Sect. ESI below we discuss the nonnegligible selection 
biases of the optical SDSS AGN sample; as a result, we 
regard the X-ray-selected AGN sample as less contam¬ 
inated by observational selection effects. Thus, we do 
not hnd convincing statistical evidence for a clustering 
dependence with L/Ledd- 

For the X-ray-selected AGNs, we still detect a 2 .Oct dif¬ 
ference between the low and high Mbh samples when we 
do not require that the subsamples have matched dis¬ 
tributions in L/Ledd- Using simple cuts in L/Ledd 
without required matched distributions in Mbh yields a 
2 . 2 ct difference in that low and medium L/Ledd AGNs 
cluster more strongly than high L/Ledd AGNs. If we 
assume that this L/Ledd AGN clustering dependence is 
real and that there is no clustering dependence on Mbh, 
this finding would contradict the observed X-ray lumi- 
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Fig. 10. — Left: HOD-derived bias parameter as a function of SMBH mass for the optically selected SDSS (blue) and the X-ray-selected 
RASS/SDSS AGN sample (red). The diamond points show a division of the full sample in Mbh 30% high, 40% medium, and 30% 
low Mbh subsamples. The boxes (with dotted error bars) indicate the combined medium and low Mbh subsamples. All data points are 
plotted at the median Mbh of the corresponding subsample. 

Right: similar to the left panel, showing the HOD-derived bias parameter as a function of Eddington ratio. For the optically selected AGN 
sample, we display the division of the full sample into 30% high, 40% medium, and 30% low (solid blue line) and the division into 40% 
high, 20% medium, and 40% low (dashed blue line). For the sake of clarity, for the X-ray-selected AGN sample we show only the division 
into 30% high, 40% medium, and 30% low, which is extremely similar to the other divisions for this sample. 


nosity clustering dependence that high Lx AGNs (and 
thus high L/Ledd) are more strongly clustered than low 
Lx AGNs. Figures H] and [S] show that the AGNs with 
the lowest L/Ledd AGN also have high Mbh- Thus, 
we conclude that differences in Mbh (in these naively 
defined subsamples) are responsible for the observed dif¬ 
ferences in the clustering of the low and high L/Ledd 
AGN subsamples. Similar results are found for the opti¬ 
cally selected AGN sample when we do simple divisions 
of the full sample without matching the distributions in 
the other parameter of interest. 

Since Mbh is estimated from the observed parameters 
FWHM and Lho, we also investigate the dependence of 
the clustering strength when we divide the full samples 
according to these parameters (see Table [T]). As before, 
we divide the sample in one parameter in such a way 
that we conserve the distribution in the other parameter 
in all subsamples. In our sample of luminous broad-line 
(X-ray-selected) RASS/SDSS AGNs, we find a 2 .Oct clus¬ 
tering difference between the lowest and highest FWHM 
samples and only a 0.6 ct difference between the lowest 
and highest Ha luminosity samples. For the optical AGN 
sample, we find a 1.9 ct difference for the low and high 
FWHM samples and a 1.5 ct difference for the low and 
high Ha luminosity samples. Both parameters are used 
to estimate Mbh, though the scaling is much stronger 
with FWHM than with Lhq (converted to L 5100 , see 
Equations [T] and [2]) ■ The similar FWHM clustering de¬ 
pendence in the X-ray and optical AGN samples provides 
further evidence that Mbh is the key parameter driving 
the observed Lx dependence of the clustering strength. 

Analysis of the two-dimensional parameter space in 
Fig. m further supports our claim that the Lx depen¬ 
dence of the clustering strength originates from an Mbh 
dependence. For the X-ray-selected sample (Fig. [8]), the 
contours of the low and high Mbh samples are extremely 
similar to the Lx contours, while the L/Ledd contours 
differ significantly from the Lx contours. We therefore 


conclude that only the black hole mass, and not a com¬ 
bination of Mbh and L/Ledd, is responsible for the ob¬ 
served weak X-ray luminosity dependence of the cluster¬ 
ing. Within the explored Mbh range, we estimate an 
average increase of the bias parameter of Ab ~ 0.7 per 
dex Mbh- This corresponds to approximately one dex 
increase in Mdmh per dex Mbh- 

If our studies had shown that the clustering strength 
correlates directly with L/Ledd but not with Mbh, then 
more strongly clustered AGNs would accrete more ma¬ 
terial. In this scenario, because highly accreting AGNs 
would lie in environments with a high density of galax¬ 
ies, galaxy mergers and galaxy-galaxy interactions would 
likely be the main triggers of AGN activity. However, we 
find that the clustering strength depends mainly on Mbh 
and not on L/Ledd- Thus, high X-ray luminosity AGNs 
do not necessarily require dense environments to accrete 
more matter. We find that, on average, AGNs in dense 
environments have higher Mbh than their counterparts 
in lower density environments. The higher Lx that are 
observed are then a direct consequence of having higher 
Mbh- 

Our findings are also consistent with a scenario in 
whic h the AGN luminos ity can vary on short timescales 
le.g.. lHickox et ani20I4f) . due to a variability in the black 
hole accretion rate (L/Ledd)- In this view, the instan¬ 
taneous luminosity of an AGN is only a weak indica¬ 
tor of the time-averaged black hole accretion rate. Since 
-^/Ledd (which reflects the amount of accreted gas) can 
rapidly fluctuate, it will not necessarily show a corre¬ 
lation with the properties of the host galaxy or DMH 
mass. On the other hand, Mbh cannot fluctuate on short 
timescales as it stays nearly constant or grows slowly. 
Therefore, the black hole mass is not subject to stochastic 
variability and correlates more tightly with host galaxy 
properties and the host DMH mass. 

In the hierarchical model of structure formation, more 
massive galaxies reside in more massive DMHs (e.g.. 
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Fig. 11.— Similar to Fig. 1101 Left: HOD-derived bias parameter as a function of absolute SDSS Mi magnitude for low (33%), medium 
(34%), and high (33%) subsamples. Right: HOD-derived bias parameter as a function of absolute rest-frame magnitude in the 5500-6800 
A band. The data points shown with a box (with uncertainties shown as dotted lines) correspond to the samples created by combining the 
low and medium absolute 5500-6800 A rest-frame magnitude subsamples. 



Fig. 12 .— Simple model prediction of the relati on between the 
dark matter halo mass and SMBH mass, using IBehroozi et al.l 
ll20il) and IConrov &: Whitd 1I2013IV This simple model does not in¬ 
clude any potential effects that an AGN can have on the host galaxy 
or host halo. The gray error bars of the correlatio n reflect the ob¬ 
servat ional errors in the data used to derive the IBehroozi et al.l 
(120131) relationship between stellar mass and M^mh- The trian¬ 
gles show our clustering results for the X-ray-selected AGN sample 
when we divide the full sample into low, medium, and high Mbh 
samples. 

IMostek et al.l 1201^ ISkibba et al.l 1201^ . More massive 
galax i es are also mor e lum inous (e.g., iTullv fc Fisheil 
Il977t iMcGaugh et al.l 1200011 . This leads to the well- 
observed luminosity d ependence of the c lustering sig¬ 
nal fo r galaxies (e.g., IZehavi et al.l 120111 : ISkibba et al.l 
120141 . If more massive galaxies also had more mas¬ 
sive bulges, the relationship between the Mbh and stel¬ 
lar ve locity dispersion of the sur r ounding galactic bulge 
le.g.. iFerrarese fc Merritd 120001 : iGebhardt et al.l 120001 : 
iKormendv fc Hoi 1201311 should then also lead to a cor¬ 
relation of Mbh and Mdmh- 
We first compare our results with a simple empirically 
motivated model, rather than a full semianalytic model. 
We take the relations bet ween galaxy ste l lar m ass and 
DMH mass at z ~ 0.3 from IBehroozi et al.l (|2013ll . Next, 
we assign each galaxy a black hole with a mass that is 
determined by the stellar mass of the host galaxy, follow¬ 


ing |Conro2EWhiti (|2013[1 (their Eq. 1). No scatter in 
the stellar mass to Mbh relationship is included in our 
treatment. 

The result of this simple model is shown in Fig. [T^ 
along with our measurements. There are discrepancies 
between the model prediction and measurements: at 
Mbh < 1O®M0 the model predicts a lower Mdmh than 
we find, though the difference is within 3cr. At higher 
MbH) the model predicts a stronger relationship between 
Mbh and Mdmh than we find. While it is understand¬ 
able how our measurements could lie about the predicted 
line, in that we measure the mass of the parent halo 
from the large-scale bias while the model predicts the 
mass of the subhalo hosting the AGN, our measurement 
at Mbh ~ IO^'^Mq falls below the prediction (a 3.7cr 
difference). 

The shallower slope seen in the data compared to the 
model here could be caused by the cumulative scatter in 
each of the correlations in the logical chain (Mdmh —t 
Msteiiar “t Mbuige —t Mbh)- We will e valu ate the predic¬ 
tions of semianalytic models in Sect. 16.41 and will show 
that they describe our findings better than the simple 
model does. 

6.2. Comparing Our Results to Other Studies 

iZhang et al.l (|2013ll and iKomiva et al.l (120131 1 study the 
projected quasar number density around galaxies. Both 
find evid ence that the cluste ring scale length cq depends 
on Mbh - IZhang et al.l (120131 1 detects this trend, however, 
at a significance level of only ^la. They also note that 
there is no clustering dependence o n optical quasa r lumi - 
nosity. Based on SDSS AGN data, iKomiva A al.l (1201311 
find a trend that tq increases with Mbh > 10® Mq. Our 
ro measurements for the Mbh subsamples (although at 
moderately different redshift ranges) agree with their re¬ 
sults when we co nsider the la uncert ainties. At lower 
black hole masses. TKomiva et al.l (1201311 do not find an rp 
dependence. This might be due to selection bia ses o f the 
AGN pop ulation, as we will sho w below in Sect. 16.41 Ad¬ 
ditionally, iKomiva et al.l (|2013ll do not find significant ro 
dependences on AGN luminosity (Asioo)- However, the 
dynamic range in their study is less than a magnitude. 
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iShen et ^ (I2009D compute the ACF of optically 
selected broad-line SDSS AGNs (DR5) at 0.4 < z < 2.5. 
They detect a ~2cr difference in the clustering strength 
when they split the AGN sample into the 10% most 
massive M^bh and the remaining 90%, in that AGNs 
with the highest Mbh are more strongly clustered. Their 
results are consistent with our findings, although their 
sample extends up to z = 2.5 and con tains more massive 
black holes at higher redshifts. When [Shen et al.l (l2009fl 
split their sample into luminous and faint quasars, 
they found that both correlation functions agree with 
each other. Only the 10% most luminous quasars show 
a larger clustering strength (at the ^2cr level) than 
the remaining 90% of the sample. This result is also 
consistent with our findings, which we will discuss in the 
next subsection (see also Fig. [TTl left). 


6.3. Why Do We Not Detect an M[ Clustering 
Dependence in the Optically Selected AGN Sample? 

Interestingly, we find very similar clustering depen¬ 
dences for the X-ray and the optically selected AGN sam¬ 
ples. However, an important question remains: if there 
are such similar dependences in the clustering strength 
with Mbh, L/Ledd, FWHM, and Ana, why do we not 
see a clustering dependence with Mi when we see a de¬ 
pendence with Lx? Within the X-ray-selected sample, 
we detect a very weak Mi clustering dependence, in that 
the faint and luminous Mi RASS AGN samples (see Ta- 
ble[T]) differ by 1.2 ct. The direct comparison between the 
Mi clustering dependence of the X-ray and optically se¬ 
lected sample in Fig. [TT] (left) shows that only the faint 
optical Mi subsamples deviate from the trend detected 
in the X-ray-selected sample. However, the difference is 
not significant considering the uncertainties of the mea¬ 
surements. 

Above z ~ 0.28 the Ha line is shifted outside the 
SDSS i-band filter. The redshift range studied here in¬ 
cludes sources above and below this redshift. Thus, be¬ 
low z ~ 0.28 the broad-line Ha line contributes flux to 
the Mi band, while above z ^ 0.28 this direct AGN lu¬ 
minosity indicator is not contributing to the Mi band 
flux. To explore if this difference in Mi contribution is 
responsible for the observed clustering difference in X- 
ray and optically selected AGN samples, we derive the 
absolute 5500-6800 A rest-frame magnitude for all X-ray 
and optically selected AGN. The band pass of this rest- 
frame magnitude is chosen in such a way that it includes 
the Ha line for all objects in our samples. We compute 
the clustering dependence as a function of this rest-frame 
absolute magnitude (Fig. [TTJ right) and find trends very 
similar to those as a function of Mi (Fig. [TTl left). Thus, 
shifting the Ha line outside the M{ band does not explain 
the difference in the clustering signals. 

Figure [13] shows that the luminous absolute magnitude 
subsamples for the X-ray and optically selected AGN 
cover a similar parameter space in the Mbh versus ab¬ 
solute 5500-6800 A rest-frame magnitude plane. The 
median values of the samples are almost identical. Thus, 
it is not very surprising that these samples have similar 
clustering strengths given their similar values of Mi, Ha, 
and absolute 5500-6800 A rest-frame magnitude. On 
the other hand, Fig. [13] demonstrates that the faint ab¬ 


solute magnitude subsamples for the X-ray and optically 
selected AGN samples are substantially different. The X- 
ray sample extends (as in all other parameters) to much 
lower magnitudes. Thus, the dynamic range of the X-ray- 
selected RASS/SDSS sample is wider than that of the 
optical AGN sample. As a consequence, the medium ab¬ 
solute magnitude subsample of the X-ray-selected AGN 
sample includes the vast majority of the optical AGN 
in the faint absolute magnitude subsample. If one com¬ 
bines the faint and medium absolute magnitude subsam¬ 
ples in the optical to roughly match the medium X-ray 
subsample, the clustering strength between both AGN 
populations agrees very well (6opt.faint-i-med. = 1.29tQ;g7; 
see Fig. HH). 

Finally, we explore the selection of the X-ray and op¬ 
tical AGNs in more detail. The RASS/SDSS sample 
is an X-ray flux-limited sample. ROSATA soft energy 
range significantly biases the sample toward broad-line 
AGNs as these sources are unabsorbed (or only mildly 
absorbed) in the X-rays. In addition to the X-ray flux 
limit, the sample is also subject to an optical flux limit 
because an SDSS spectrum is required for the charac¬ 
terization of a source. However, the X-ray and opti- 
cal luminosities of broad-line AGNs correlate well (e.g., 
Tananbaum et al.llI979l: iGrupe et al.ll20ldl: iLusso el's] 
2010(1 . Since ROSAT has limited sensitivity, it is able to 
detect only the X-ray-brightest AGNs. These objects are 
typically also bright in the optical. Thus, the impact of 
the optical flux limit on the sample selection is moderate. 

The optical SDSS AGN selection, on the other hand, 
is much more complex. The target selection relies on 
several different selection techniques. The main criterion 
is a color cut and an optical flux limit. Broad-line AGNs 
with a certain amount of extinction will not meet the 
color-cut criterion. Additional objects are included in the 
sample due to their radio emission. Most importantly, an 
absolute magnitude cut is applied. 

Figure [14] (left) shows that the optically selected AGN 
sample is not a flux-limited sample, although the sample 
selection considers only objects with 15.0 < i < 19.1. 
The absolute Mi cut (applied to prevent substantial con¬ 
tamination by host galaxy light) causes a serious defi¬ 
ciency of objects above the lower i-band limit off < 19.1. 
The X-ray-selected AGN sample lFig.[T4l right) includes 
objects down to almost the same apparent i-band mag¬ 
nitude, independent of redshift. Various SDSS papers 
discuss how the optical SPS S AGN sample is not a uni¬ 
formly selected sample le.g.. [Richards et ^l2006bll . We 
conclude that the complex optical AGN selection might 
be the origin of the clustering differences between the X- 
ray and optically selected AGN samples. The strongest 
differences are detected between the samples of the low¬ 
est median redshift; this is also the range in which the 
optical AGN sample deviates most from a clean, flux- 
limited sample. We therefore consider the X-ray-selected 
RASS/SDSS AGN sample as the more complete sample 
and use only their clustering properties in the next sec¬ 
tion. 


6.4. Comparison to cosmological simulations 

Hydrodynamica. 1 __ iBrvan et al.l 12014 

iHirschmaim et ah! 12014 iS teinbor n et alJ 2015(1 and 
semianalytical (e.g.. iMalbon et al.l l2007( : iMarulli et al.l 
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Fig. 13.— Black hole mass for optically selected SDSS (left) and X-ray-selected RASS/SDSS (right) AGNs vs. absolute 5500-6800 A 
rest-frame magnitude. Different colors denote the division of the full samples into faint (33%), medium (34%), and luminous (33%) objects 
with respect to the rest-frame absolute magnitude. The black diamonds show the median absolute magnitude and median Mbh of fh® 
subsamples. We also show linear regression line fits to the full samples as dotted lines. 




Fig. 14.— Apparent SDSS f-band magnitude vs. redshift for optically selected AGNs (left panel) and X-ray-selected AGNs (right panel). 
The horizontal dashed lines represent the lower and upper f-band limits for the selection of optical SDSS AGNs. 


I2008I iFanidakis et al.l 120111 : iBensonI l2012t) state-of- 
the-art simulations have reached a level of complexity 
such that they can well predict observed quantities 
such as luminosity and stellar mass functions (e.g., 
HopUns et al.l[200a iDeeraf et al.ll2010l IFanidakis A all 
201211 . The observed clustering properties of galaxies 
and AGNs provid e additional impor tant constraints for 
such models le.g.. iBonoli et al.ll2009fl . 

Theoretical models of galaxy formation that include 
prescriptions for modeling BH growth can provide impor¬ 
tant constraints on the fueling modes of AGN. For exam¬ 
ple, models in which AGN activity is triggered by major 
mergers between galaxies or secular processes within the 
host galaxy (e.g., disk instabilities) predict a well-defined 
correlation function for luminous quasars, with an aver¬ 
age DMH niass of Mpmh ^ 10^^ h~^ Mq (|Bonoli et al.l 
120091120101 : IFanidakis et al.l[2013H . However, the inclu¬ 


sion of an additional mode of grow th that is linked to dif- 
fuse gas accretion in massive halos (|Fanidakis et al.lF2012[) 
predicts a much higher average halo mass for moderate- 
luminosit y AGNs, in contrast to their luminous coun¬ 
terparts (|Fanidakis et al.l[2013D . Gomparing this specific 
picture of AGN clustering with the available observations 
could potentially provide important insights into the cor¬ 
rect modeling of the growth of black holes in galaxy- 
formation models. 

Cosmological simulations and their predicted cluster¬ 
ing properties can be also be used to test the impact of 
sample selection biases for observed samples. As an ex¬ 
ample, let us assume that in an unbiased sample there 
is no clustering dependence as a function of a certain 
AGN parameter. A particular observational selection 
bias might cause a clustering dependence to be measured. 
As selection biases can be included in the cosmological 
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Fig. 15.— Eddington ratio vs. mass of the SMBH of simulated 
(color) and observed (gray) X-ray-selected AGNs. The color de¬ 
notes the number of AGNs above the RASS flux limit divided by 
the total number of AGNs in the simulation, after applying the 
RASS flux limit. The gray data points indicate the location of the 
observed X-ray-selected RASS/SDSS AGN sample (see Fig.[5jl. 

simulation, one can therefore compare the predictions for 
the unbiased and observed (including selection biases) 
AGN samples. Such knowledge is crucial to avoid wrong 
interpretations about the underlying physics when using 
clustering studies. Thus cosmological simulations are not 
only important to test our theoretical understanding of 
the observed universe but also to explore the impact of 
observational sample biases. 

Here w e use the semiana lytical galaxy formation model 
GALFDRM (iCole et al.ll2000fl to compare with our observed 
AGN clustering signals. iFanidakis et al.l (|201lL I2012L 
[Mil modeled SMBH in these simulations, using an 
identical model in all three papers. Only the cosmol¬ 
ogy has chan ged from WMAPl (iii the first two papers) 
to WMAP7 (IFanidakis et al.l|2013ll . In IFanidakis h all 
(|2013ll . the authors incorporate two modes of AGN ac¬ 
cretion in the GALFORM model. The first is the starburst 
mode (cold accretion), in which accretion onto the BH 
is tightly coupled to the mass of the cold gas available 
in the galaxy, which also contributes to star formation. 
This mode occurs when the host galaxy encounters a ma¬ 
jor galaxy merger, a minor merger, or a disk instability. 
These phenomena can take place in a wide range of halo 
masses. However, observable AGN luminosities within 
the range defined in this study are ty pically produced 
in DM Hs with Mdmh < h~^ Mq (IFanidakis et al.l 

IMl). 

The second is the hot-halo mode (hot accretion), which 
occurs in DMHs with higher masses. A fraction of the 
AGN energy output is used to heat the gas in the DMH 
and suppress cooling of the gas. The SMBH accretes 
directly from the diffuse gas in the DMH. This mode 
regulates the black hole accretion in very massive ha¬ 
los with Mdmh ^ 10^^/i“^Mq. For more details on 
both modes and their modeling, we refer the reader to 
IFanidakis et al.l (|2013ll and the references within. 


The model of IFanidakis et all (I2013D is successful in 
reproducing several observational constraints, including 
the observed luminosity functions of galaxies and AGNs 
over a wide range of redshifts (jFanidakis et al.l 120111 
l2012f) . In addition, they can reproduce reasonably well 
the observed MBH-Mcai.buige relationship and the global 
and active black hole mass functions. 

The primary motivation for running the simulation is 
to explore how AGN clustering depends on the mode of 
AGN accretion, as described in the two modes mentioned 
above. Here we can thus compare our observational re¬ 
sults directly to theoretical predictions. In addition, we 
can explore the effect of the RASS flux limit selection on 
o ur observational s ample . 

IFanidakis et al.l (l2()13f) present the typical DMH 
masses for AGNs observed in a redshift range of z = 0- 
1.3. Here we repeat this simulation in a redshift range 
of z = 0.16-0.36, identical to the redshift range used for 
the observed X-ray and optically selected AGN samples. 
For each galaxy, the simulation allows one to identify 
when the AGN is active. An object is included in the 
AGN sample if its central engine is active in at least one 
of the five logarithmically spaced snapshots output be¬ 
tween z = 0.16-0.36. AGN properties such as Mbh, Aboi, 
T/Tedd, and redshift are available for each of the differ¬ 
ent snapshots. The parameter Lboi is calculated assum¬ 
ing that the accretion flow forms a ge ometrically thin disk 
for re latively high accretion rates (jShakura fc Sunv;r^ 
Il973f) . At lower accretion rates, a geometrically thick 
disk or an a dvection-dominated accretion flow is mod¬ 
eled (ADAF. [Naravan fc Yilll993l . 

X-ray luminosities are calculated directl y from £bni by 
applyi ng the bolometric correction from [Marconi et al.l 
(12004) . We consider an object to be an AGN if its intrin¬ 
sic 2-10 keV rest-frame luminosity is L 2 -iokeV > lO'*^'® 
erg s“^. The sample contains X-ray unabsorbed and 
absorbed AGNs. Ba sed on th e lumi nosity, we use the 
empirical formula of iHasingeil (j2008l l to determine the 
probability that the object is absorbed. We then clas¬ 
sify statistically if an object is unabsorbed or absorbed. 
Below, we will refer to this AGN sample as the “all sim¬ 
ulated AGN” sample. More recent estimates of the ob¬ 
scured fraction of AGNs as a function o f luminosity have 
been presented by iMerloni et al.l ()2ni4 l , who show that 
the probability that an AGN at high l uminosities is ab - 
sorbed is higher than that suggested bv lHasing^ (1200811 . 
We test the impact of a higher obscured fraction at high 
luminosities and find insignificant differences as far as 
clustering measurements are concerned. We therefore re - 
tain the empiri cal obscuration relation of lHasing^ (j2008f l 
as presented in IFanidakis et ahl (j2013f l. 

Since RASS is essentially a flux-limited survey (see also 
Fig.im right), only AGNs above a certain X-ray flux are 
contained in the RASS AGN catalog. However, due to 
the characteristics of the RASS, different regions of the 
sky have different exposure times. To mimic the RASS 
flux selection, we also apply an average flux cut to the 
simulated AGN sample by using the X-ray luminosities 
and redshifts of objects in the simulation. We refer to 
this sample as the “RASS-selected AGN sample.” Its 
location in the Mbh versus A/Tedd plane is shown in 
Fig. dSl (colored grid points). This sample also includes 
the effect that RASS is sensitive only to unabsorbed X- 
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Fig. 16.— Left: comparison between the predicted (lines) and observed (data points) RASS AGN typical host dark matter halo mass 
as a function of hard (2-10 keV) X-ray luminosity. The red line shows the prediction of a simulated RASS flux-limited sample (and log 
(L/Lbdd) > ~2), and the black dash-dotted line shows the expected dark matter halo mass for the total (unbiased) AGN sample with 
^2— lOkeV 10“^^ ® erg s"*^. We recalculate the relations after adding a random error from a Gaussian distribution with cr = 0.3 dex to the 
simulated SMBH masses. We run 50 such realizations; a random subset of these is shown here as orange lines. The data points are plotted 
at the c orre sponding median T2-l0keV values of the observed low and high AGN subsamples (see Table[T] using the conversion presented 
in Sect. [21^ . Middle: similar to the left panel, here showing the predicted and observed halo mass as a function of Mbh- Right: similar 
to the left panel, here showing the predicted and observed halo mass as a function of Eddington ratio (T/Tedd). 


ray AGNs. However, since the luminosities of the RASS 
AGNs are high, the simulation classifies the vast majority 
of the objects as X-ray unabsorbed. 

In the same plane, we also show the observed X-ray- 
selected broad-line RASS/SDSS AGN sample. In gen¬ 
eral, the simulation agrees well with the observations, in 
that sources lie in a similar region of this space. How¬ 
ever, the simulation slightly overestimates the number of 
AGNs with Mbh > 10®-^ Mq and log (A/Aedd) < —2. 
Such objects should be detectable given the RASS flux 
limit, but they do not exist in large numbers (at the ob¬ 
served redshift range). However, these are the objects in 
the simulation for which the accretion flow is modeled 
with an ADAF. These AGNs are observationally associ¬ 
ated with rad io-bright, mechanical feedback-dominated 
SMBHs (e.g.,[Ghurazov et al.l 120051: iHickox et al.l 1200^ 
iSmolcid 120091) . The optical spectra of these objects do 
not c ontain any high excitat ion lines typical of AGNs 
(e.s.. [Best fc Heckmanl [20I2D . Thus, one might argue 
that these objects would not be classified as AGNs in an 
SDSS spectrum. We therefore limit the simulated RASS- 
selected AGN sample (and the predictions in the follow¬ 
ing paragraphs) to those sources with log (A/Aedd) > 
- 2 . 

Despite these minor differences, we calculate the bias 
parameter for the simulated AGNs and the observed 
RASS-selected AGN sample. The bias is calculated in 
bins of 0.2 in l og (MBn/Mr?^) and log (A/Aedd) using 
Eq. 4 given in iFanidakis et al.l (I201.1D (see also Eq. 8 
in this paper). We follow the same procedure as for 
the observed AGN samples. Thus, the bias values are 
converted into typical DMH masses (log Mn^^^^). Since 
the sim ulation uses the conversion given by iSheth et al.l 
(1200 ID without applyi ng the improved fitting formula 
bv [Tinker et al.l (j2005D . we recompute for the observed 
Ax, Mbh, and A/Aedd AGN subsamples in Table [T] 
Mn/,P„ base d only on the ellipsoidal collapse model of 
ISheth et al.l (|200lD . This allows us to directly com¬ 


pare the predictions with the observations, as shown in 

Fig. HH 

In the case of the simulated RASS-selected AGN sam¬ 
ple, we also assign to every BH mass a random error from 
a Gaussian distribution with a = 0.3 dex. We recalcu¬ 
late all of the relations (red lines) shown in Fig.[Tn]for 50 
realizations and show a subset of these as orange lines. 
This exercise does not only allows for a more realistic 
comparison to the data, but also allows us to explore the 
impact of the intrinsic scatter of approximately 0.3 dex 
in the Mbh estimate of the observed sample. Fig. [11] 
shows that the 0.3 dex uncertainty in the BH mass es¬ 
timates has no significant impact on the prediction of 
the correlations shown. Only at very low Mbh and high 
A/Aedd does the scatter between different realizations 
increase. However, this is caused by the low number of 
objects with such properties in the simulations. 

For the all-simulated AGN sample (black line), the sim¬ 
ulation predicts a positive clustering dependence on X- 
ray luminosity only above log {L 2 - io/[ergs ^]) - 44. 
This might explain why studies of moderate X-ray lumi¬ 
nosity do not see a clustering dependence on luminos¬ 
ity and conflicting results are presented in the literature. 
The RASS selection results in a luminosity dependence 
that extends down to log (A 2 _io/[ergs“^]) ^ 43 (red 
line). The data and the predictions for the RASS sample 
agree remarkably well. 

For the Mbh clustering dependence (Fig. [TH middle), 
the simulation predicts a strong clustering dependence 
with Mbh for the all simulated AGN sample; a steady in¬ 
crease of typical DMH mass with Mbh is found (see also 
Fig. 7 in IFanidakis et al1l20I2D . Above Mbh ~ 10® Mq 
the correlation is even stronger. The RASS flux limit se¬ 
lection of the all-simulated AGN sample has a moderate 
impact on the Mbh clustering dependence. It amplifies 
the correlation at Mbh ^ 10® Mq and weakens it below 
AIbh 10® Mq compared to an unbiased AGN sam¬ 
ple. The model shows that, in flux-limited samples be- 
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low Mbh ^ 10® M 0 , no clustering correlation with Mbh 
should be detectable, although an unbiased AGN sample 
will still show a very weak Mbh clustering dependence in 
thi s Mbh range. This is consistent with the result found 
bv iKomiva et al.l (I2013I1 . 

The observed Mbh clustering dependence agrees rea¬ 
sonably well with the model prediction. The data hint 
that the observed clustering dependence might not be 
as strong as predicted by the simulated RASS-selected 
AGN sample. 

The predicted T/Tedd clustering dependence also 
matches the observations well (Fig. (THl right). Below 
L/Ledd ~ 10“^ the model predicts higher DMH masses 
than are currently detected for AGNs, but our data do 
not probe Eddington ratios that low. The model predicts 
that, at low L/Ledd values, the RASS flux limit selection 
should cause a very strong negative T/Aedd dependence 
to the clustering amplitude. The data point with the 
lowest A/Aedd in the RASS/SDSS AGN sample does 
not allow us to verify or reject this prediction. Although 
our samples contain objects down to A/Aedd ~ 10“^ 
(see Fig. SI), AGN samples with even lower A/Aedd are 
needed to critically test the prediction of high Mdmh at 
low A/Aedd- 

Such tests between observations and state-of-the-art 
simulations offer a unique opportunity to constrain the 
physical mechanisms included in galaxy and AGN forma¬ 
tion models. For example, a mechanism for the hot-halo 
mode in which black holes with low A/Aedd and high 
Mbh reproduce the real RASS-selected sample might 
agree better with the observed AGN clustering depen¬ 
dences. Thus, a channel has to be found to remove 
the simulated objects with Mbh ^ 10®'® Mq and log 
(A/Aedd) ^ —2. In the current model, the hot-halo 
AGNs arise from the necessity for AGN feedback in mas¬ 
sive DMHs. The accretion rate in this mode is calcu¬ 
lated indirectly from the cooling properties of the host 
halo and corresponds to an accretion flow that is able 
to reproduce powerful AGN outflows that suppress gas 
cooling and star formation. The tuning of the free pa¬ 
rameters in this calculation is done by fitting the model 
to the observed galaxy luminosity function at 2 : = 0. 
Obviously, more observational data are needed to con- 
train the model. The comparison between the simulated 
and real RASS-selected AGN samples in Figs. [T5] and ITbl 
could be used as additional constraints for constraining 
the free parameters in the model. 

When cosmological simulations such as the one 
presented here are adjusted to match the observed AGN 
clustering properties, to be considered successful they 
still have to match other observational constraints, such 
as the luminosity function of galaxies and AGNs at 
different redshifts. Changes in the physical treatment of 
AGN accretion might also solve other challenges such as 
the observed deficiency of low-luminosity AGN at high 
redsh ifts when compared to predictions from simulations 
(e.g., iMivaii et al.ll201,^ . Including various constraints 
from galaxy and AGN clustering measurements has thus 
the power not only to improve the simulations but also 
to enhance our general understanding of AGN physics 
and AGN and galaxy coevolution. 

7. CONCLUSIONS 


Motivated by the detection of a weak X-ray luminosity 
dependence of the clustering strength of AGNs in the 
first paper of this series, here we explore the physical 
origin of this dependence. Using the optical spectra of 
our soft X-ray-selected (RASS/SDSS) luminous, broad- 
line AGN sample at 0.16 < 0.36, we estimate black hole 
masses and Eddington ratios and calculate the clustering 
dependence on each. 

Since Mbh and A/Aedd are correlated, we create sub¬ 
samples in Mbh and A/Aedd that have matched dis¬ 
tributions in the other parameter of interest. We com¬ 
pute the clustering strength for the subsamples and find a 
weak clustering dependence with Mbh and no significant 
correlation with A/Aedd- Various adjustments in how 
the subsamples are created do not change the results. 
We also study the clustering of the observed parameters 
(luminosity and FWHM of the broad Ha line) that are 
used to derive Mbh- We find a weak FWHMho cluster¬ 
ing correlation in that AGNs with low FWHMna are less 
clustered than their high FWHMna counterparts. 

We also study the clustering properties of an optically 
selected SDSS AGN sample. This sample has 29% of its 
objects in common with the X-ray-selected RASS/SDSS 
AGN sample. We detect the same trends as found for the 
X-ray AGN sample with respect to Mbh, A/Aedd, Ahq,, 
and FWHMho- The X-ray and optically selected AGN 
samples show divergent clustering signals for the lowest 
Mi (absolute Aband magnitude) samples. We argue that 
this is caused by various complex selection effects in the 
optical sample and that the X-ray-selected RASS/SDSS 
AGN sample is more uniformly selected. 

From our correlation function measurements, we con¬ 
clude that Mbh is the origin of the observed weak X- 
ray luminosity clustering dependence. The confidence 
contours of our AGN HOD modeling parameters further 
support this finding, as low and high Mbh and Ax sam¬ 
ples show extremely similar contours, while the low and 
high A/Aedd sample contours differ significantly from 
the Ax contours. 

The Mbh clustering dependence is detected at a sig¬ 
nificance level of 2.7cr. In both the X-ray and optical 
AGN samples, the highest clustering strength is found 
for AGNs with the 30% highest Mbh- Thus, at a red- 
shift range of 0.16 < z < 0.36 luminous broad-line AGNs 
with more massive Mbh reside, on average, in more mas¬ 
sive DMHs. In this context, the DMH mass refers to the 
single largest (parent) halo mass. 

Since the observed clustering strength does not depend 
on A/Aedd in a statistically significant way, AGNs with 
high accretion rates do not require large-scale dense en¬ 
vironments with high galaxy density. This provides ev¬ 
idence that major or minor mergers play only a limited 
role in the AGN accretion processes in the low-redshift 
Universe. Internal processes such as disk instabilities 
could be the dominant AGN triggering mechanism at 
late cosmic times. 

Empirically motivated models that include simple 
monotonic relationships between Mdmh, -^stellar, and 
Mbh and without AGN feedback or treatments of AGN 
triggering do not match well the observed relationships 
between Mbh and Mdmh derived from clustering. Thus, 
we use a semianalytical model to study the predicted 
clustering dependences as a function of Ax, Mbh, and 
A/Aedd- This model included AGN feedback and two 
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modes of AGN accretion and fits the observed Mbh vs. 
AIdmh relation better than the simple model. Com¬ 
paring the simulated full AGN sample with a simulated 
RASS flux-limited AGN sample, we show that observa¬ 
tional selection effects moderately change the expected 
clustering dependences. The simulation predicts that 
higher Mbh are found in more massive halos, and the 
predicted correlations with Lx, Mbh, and L/Ledd agree 
reasonably well with our observations presented here. 

Clustering measurements with small uncertainties offer 
a unique opportunity, after considering the effects of se¬ 
lection biases, to identify missing physics in our current 
understanding of the galaxy and AGN evolution by com¬ 
paring theoretical predictions with observations. Thus, 
future clustering measurements as a function of various 
galaxy and AGN parameters will deliver needed addi¬ 
tional constraints that will improve upon and distinguish 
between current theoretical models. 
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